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Preface 



Logistics of transport systems is a key driver for the growth of whatever economy. 

Freight transport allows production systems or common citizens to receive or send 
materials or finished goods required by processes as well as by everyday life 
activities. Passenger transport, both public and private, allow people saving time 
for their transfers and ensure high level of mobility. Materials and people journeys 
fulfill economy and society expectations. 

However, the overall transport sector accounts worldwide for more than half of 
global liquid fossil fuels consumptions which, in turn, is responsible for a nearly 
quarter of the world’s energy-related CCL emissions, more than 80 % of air 
pollution in the cities and about 1.3 million of fatal traffic accidents per year. 
Negative effects represent ‘external costs’ paid by unaware societies and modern 
economies. 

Costs of externalities account worldwide for more than 10 % of the GDP with an 
increasing trend. The European Environment Agency (EEA) and the United 
Nations Environment Programme (UNEP) defined the ‘Avoid-Shift-lmprove’ 
(ASI) strategy to tackle the increasing of externalities while EU Commission 
(Directorate General for Mobility and Transport) established in 2011 a roadmap 
that will lead to the internalization of external costs within 2020. Research 
programs and strategic actions on sustainable development of smart cities are 
focusing on smart mobility of goods and citizens due to the relevant 
environmental, social and economic costs of logistics. 

Internalization of cost of externalities gives rise to new logistics cost estimates and 
functions which managers, researchers, lecturers and students should refer in 
facing with logistics issues. Under this purpose the present book has been 
conceived. 

The book focuses on freight transports of industrial production systems. The most 
used keywords are as follows: sustainable logistics, freight transport, 
internalization of external costs, environmental cost, social cost, inventory 
management. Economic Order Quantity — EOQ, logistics cost function, loss factor 
of transport, Sustainable Order Quantity — SOQ, transport means selection, 
stochastic variability of product demand, stochastic variability of supply lead time, 
sensitivity analysis, finished vehicle logistics, inland waterways, automotive 
supply chain, spare parts, repair policy. 

The book has been subdivided into three main parts, organized as introduced 
below. 




vi Preface 



Chapter 1 provides a taxonomy of external cost figures as well as data set enabling 
the reader to perform reliable estimates of freight transport external costs. To this 
purpose, a full scale case study is developed. 

Chapter 2 describes a new sustainable inventory management model whose cost 
functions include externalities. The classical ‘Economic Order Quantity’ model is 
re-formulated and the new concept of Sustainable Order Quantity (SOQ) is 
defined. 

Finally, in Chap. 3 the SOQ model is formulated for different inventory 
management applications referred to both deterministic and stochastic production 
environments. Numerical examples are provided. 

We would like to thank our colleagues, both academics and professionals from 
service companies, our students, and the editors at Springer for their valuable and 
helpful support. 
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1 Internalization of External Costs of Freight 
Transport 



Abstract 

Accidents, global warming, congestion, air pollution and noise are examples of 
negative effects related to the transport activities that generate costs not fully 
borne by the transport users and hence not taken into account when they make a 
transport decision: these are the so called external costs. The internalization of the 
external costs of transport has been an important issue for transport research and 
policy development for many years worldwide. This Chapter, starting from an 
overview of the transport sector statistics and of transport external costs 
internalization in Europe, gives a taxonomy of the main transport external costs; 
moreover, the state of art on the cost estimation methodologies is briefly 
introduced. A case study from the finished vehicle logistics in the automotive 
sector is presented. Results show the potential external costs reduction due to the 
better environmental and social performance assured by the modal shift from 
road toward inland waterways transport. 



Keywords : Sustainable logistics, Freight transport, internalization of external 
costs, Environmental cost, Social cost, Finished vehicle logistics, Inland 
waterways 

1.1 Overview on the Transport System and the Legislative 
Context 

The transport sector, including the movement of people and goods by cars, trucks, 
trains, ships, airplanes, and other vehicles, is a key drive for the European Union 
Countries economic growth. It accounts for about the 5 % of the EU28 Gross 
Value Added (GVA) and employs about the 5 % of the total workforce in the 
EU28 [1]. In 2012, freight transport activities amounted to 3768 billion [t-km] 
while passenger transport ones to 6391 billion [p-km]. Figure 1.1 shows the 1 995— 
2012 data of the Gross Domestic Product (GDP), of the freight transport and of the 
passenger transport in EU28 [1] (year 1995 values: 8012 billion [€] for the GDP, 
3.07 billion [t-km] for the freight transport and 5.37 billion [p-km] for the 
passenger transport). 
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Fig. 1.1 Gross Domestic Product, passenger and freight transport trend from 1995 to 2012 in 

EU28 [1] 



Unfortunately transport sector, characterized mainly by fossil fuel driven motor 
vehicles, gives rise to negative effects. 

This sector, worldwide, is responsible of: 

• more than half of global liquid fossil fuels consumption; 

• nearly a quarter of the world’s energy-related CCF; 

• more than 80 % of the air pollution in cities in developing countries; 

• more than 1.27 million fatal traffic accidents per year; 

• chronic traffic congestion in many of the world’s urban areas. 

These negative effects cause costs, which can add up to more than 10 % of a 
country’s Gross Domestic Product, paid by the society and are likely to grow, 
primarily because of the expected growth of the global vehicle fleet. The 
continuation on a business-as-usual path will result in an increase of the global 
vehicle fleet from around 800 million to between 2 and 3 billion by 2050. Most of 
this growth will be concentrated in the developing countries. Furthermore, it is 
also expected an exponentially growth for the aviation sector (mainly due to the 
developing countries) and a growth by up to 250 % of the carbon emissions from 
shipping [2], 

The shift toward a green transport is needed. EEA and UNEP proposed an 
holistic strategy, called the Avoid-Shift-lmprove (ASI) strategy [2], to reach this 
goal. The ASI strategy aims at: 

1. avoiding or reducing the number of journeys taken; 

2. shifting to more environmentally efficient forms of transport; 

3. improving vehicle and fuel technology. 
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In this context it can be included the external costs internalization strategy. 

The European Commission focused the attention on the external costs of the 
transport for many years and, in 1995, defined the transport externalities as 
follows [3]: 

“Transport externalities refer to a situation in which a transport user either does 
not pay for the full costs (e.g. including the environmental, congestion or accident 
costs) of his/her transport activity or does not receive the full benefits from it.” 

The aim of the external costs internalization is the integration of these costs in 
the decision making process of the transport users [4] : 

• directly: through, for example, command and control measures; 

• indirectly: through market-based instruments providing the right 
incentives to the transport users such as, for example, taxes, charges and 
emission trading; 

• by combinations of these basic types: for example, existing taxes and 
charges may be differentiated by the EURO emission classes of vehicles. 

The use of market-based instruments is generally regarded as the best strategy 
to limit the negative effects of the transport requiring, however, a detailed and 
reliable estimation of external costs. In order to better define the external costs it is 
important to highlight the difference between: 

• social costs: reflect all costs occurring due to the provision and use of 
transport infrastructure (i.e. wear and tear costs of infrastructure, capital 
costs, congestion costs, accident costs, environmental costs); 

• private (or internal costs): directly paid by the transport users (i.e. wear 
and tear and energy cost of vehicle use, own time loss costs, transport 
fares, and transport taxes and charges). 

As aforementioned, the European Commission has pointed out the objective to 
charge the vehicles for the external costs they generate since 1995 [3]. European 
Directive 1999/62/EC [5] (also called Eurovignette-Directive) and its amendment 
[6] are consistent with this goal. European Directive 1999/62/EC did not includ 
all the transport means but was limited to vehicle taxes, tolls and user charges 
imposed on heavy duty vehicles (HD Vs) aiming at the harmonization of levy 
systems and at the establishment of a fair mechanism for charging the 
infrastructure costs on vehicles using them. The spatial scope of this directive was 
the Trans-European Transport Networks (TEN-T), a planned set of road, rail, air 
and water transport networks to improve the transport sector performance in the 
European Union. The Directive already recognized, in a general way, the 
possibility to address a certain amount of the toll revenues to environmental 
protection activities but the main destination of the tolls revenues was only the 
recovering of the infrastructure costs (costs of construction, operation and 
maintenance). By adopting only the user pays principle, this Directive failed in 
recognizing also the polluter pays principle: all road users were considered alike 
without considering for example the different congestion or pollution they caused. 
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Furthermore, the spatial limitation to the TEN-T may cause a traffic shift towards 
not charged networks. In 2011, the Council adopted the new Eurovignette- 
Directive [6] acting on all Member States’ motorways and not only on the TEN-T. 
Each Member State may define tolls composed of an infrastructure charge that 
considers also the negative effect of traffic congestion, and/or an external-cost 
charge related to traffic-based externalities (e.g.: air and noise pollutions). Only 
suggestions and not obligations are provided regarding the use of the revenue from 
infrastructure and external costs. 

The internalization of the external cost is treated also in the EU White Paper in 
2011 [7]: this document comprises 40 initiatives to be actuated within 2020 in the 
EU. The ‘smart pricing and taxation’ initiative is divided into two phases. The first 
phase, up to 2016, expects to start with a mandatory infrastructure charging for 
HD Vs and to proceed with the internalization of the external costs for all modes of 
transport. The second phase, from 2016 until 2020, expects to implement a full 
and mandatory external costs internalization for road and rail transport and to 
examine a mandatory internalization of the external costs on all European inland 
waterways network. The mandatory external costs internalization could also help 
achieving other objectives included in the White Paper such as shifting: (i) 30 % 
of road freight over 300 [km] to other modes such as rail or waterborne transport 
by 2030 and more than 50 % by 2050; (ii) the majority of medium-distance 
passenger transport from road to rail by 2050. 



1.2 A Taxonomy of External Costs 

According to the most recent estimates, the total external costs of transport in the 
EU27 countries (with the exception of Malta and Cyprus but including Norway 
and Switzerland) in 2008 have been estimated at about 500 billion [€], excluding 
congestion, and at about 700 billion [€] including congestion. The GDP in EU27 
in 2008 was about 12.5 quadrillion [€]: the total impact of externalities amounted 
to 5-6 % of GDP. Fig. 1.2a shows that accidents, congestion, climate change and 
air pollution represent 86% of total costs. Moreover (see Fig. 1.2b), the road 
sector generate 93 % of total external costs, rail accounts for 2 %, the aviation 
passenger sector 4 % (only continental flights), and inland waterways 0.3 % [8]. 

In the following, for each cost category, the type of cost figures considered and the 
methodologies for their estimation are pointed out. 



1.2 A Taxonomy of External Costs 5 





Fig. 1.2 External costs in EU27 in 2008 [8] 



Accident Costs 

The external accident costs are related to the costs not covered by the insurance 
premium such as, for example, pain and suffering caused by the traffic accidents. 
The best approach to estimate the marginal accident external costs is the bottom- 
up methodology. The main assumption of this approach is that when the driver of 
an additional vehicle joins the traffic exposes himself/herself to the average 
accident risk. This average accident risk can be estimated through the statistical 
relationship between the number of accidents involving a given vehicle class and 
the traffic flow observed in the previous years. The costs related to the accident 
risk are: 

• the expected costs, for the person exposed to the risk, of death and injury 
due to an accident; 

• the expected costs for the relatives and friends of the person exposed to 
the risk; 
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• accident cost for the rest of the society (material costs, police and 
medical costs, output costs). 

The concept of the Willingness To Pay (WTP) for safety is used to evaluate the 
first two cost elements focusing on the Value of a Statistical Life (VSL) [4]. The 
estimates on the VSL generally come from studies where participants to these 
studies quantify own WTP for the reduction of the accident risk. These estimates 
are different across countries, age groups and also differ from the risk analyzed: in 
fact, the expected number of life years lost differs among different risk cases. 

Several approaches could be adopted to quantify the share of the external costs 
in total accident costs taking into account what is already covered by the insurance 
of the person exposed to risk [4]. 



Climate Change Costs 

Climate Change (or Global Warming) impacts of transport are mainly related to 
the emissions of the greenhouse gases such as carbon dioxide (CO 2 ), methane 
(CH 4 ) and nitrous oxide (N 2 0). In the case of aviation, at high altitude, also other 
emissions (water vapor, sulfate, soot aerosols and nitrous oxides) have an impact 
on global warming. Several methodologies are available to estimate the climate 
change costs for the different transport modality: the state-of-art approach for 
evaluating this externality is the damage costs approach called Impact Pathway 
Approach (IPA) characterized by the following main steps: 

• quantification of the GHG emission factors, in [(C02) cq ] for different 
vehicles; 

• valuation of climate change costs per tonne of [(C0 2 ) cq ]; 

• calculation of the marginal climate change costs for different vehicle and 
fuel types. 

The damage cost approach and the abatement cost approach are the two main 
methodologies evaluating the cost of the GHG emissions [4]. The first one 
evaluates the total costs supposing that no efforts are taken to reduce the GHG 
emissions; the second one evaluates the costs of achieving a certain amount of 
emissions reduction. Between the two methodologies the abatement cost approach 
is preferred, although the damage cost approach is desirable from a scientific point 
of view; at the same time it is characterized by a high uncertainty mainly because 
it is not possible to identify and to evaluate many risks related to future climate 
change costs. The abatement costs approach mainly reflects the willingness-to 
pay, of a society, for a certain abatement level of the emissions. In the abatement 
costs evaluation, usually two different targets are considered in Europe: 

1. EU Greenhouse gases emissions reduction target for 2020 (corresponding 
to a cut of 20 % of GHG emissions compared to 1990 levels, “low 
scenario”); 
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2. a longer term target for keeping concentration of C0 2eq in the atmosphere 
below 450 [ppm] (thus keeping global temperature rise below 2 [°C] 
relative to pre-industrial levels [9], “high scenario”). 

Congestion Costs 

The concept of congestion externalities is easily understandable but difficult to 
quantify. A road network user affects, by his/her decision to use a certain network 
for driving between two different destinations, the utility of all other users who 
want to use the same network. The utility loss, aggregated over all those other 
users, is the negative external effect of the respective user’s decision to move 
between the same destinations. The utility loss is translated into costs considering 
the willingness to pay for avoiding this utility loss. Thus, the external effect is 
measured in terms of a monetary amount per trip. 

The update of the unit values for congestion costs, suggested by the last 
“Handbook on the external costs” commissioned by the European Commission 
[4], is based on the aggregated approach of the FORGE model used in the 
National Transport Model of the United Kingdom [10]. 

Air Pollution Costs 

Air pollution costs are mainly due to the emission of air pollutants such as 
particulate matter (PMjo, PM 25 ), nitrogen oxides (NO x ), sulfur dioxide (S0 2 ), 
ozone (0 3 ) and Volatile Organic Compounds (VOC). The following effects are 
related to this externality: 

• health costs. Impacts on human health due to the aspiration of fine 
particles (PM 2 5 /PMi 0 , other air pollutants). In addition, also Ozone 
(0 3 ) has impacts on human health; 

• building and material damages. Mainly two effects have the most 
impact: (1) soiling of building surfaces/facades mainly through 
particles and dust; (2) degradation on facades and materials through 
corrosive processes due to acid air pollutants like NO x and S0 2 ; 

• crop losses in agriculture and impacts on the biosphere. Acid 
deposition, ozone exposition and S0 2 damage crops as well as forests 
and other ecosystems; 

• costs for further damages for the ecosystem. Eutrophication and 
acidification due to the deposition of nitrogen oxides as well as 
contamination with heavy metals (from tire wear and tear) impact on 
soil and groundwater. 

The unit cost estimation of the air pollution for the different transport 
modalities follows the already mentioned Impact Pathway Approach aiming, in 
this case, at the quantification of the impact of the emissions on human health (the 
major effects), environment, economic activity, etc. [4], The key steps of the IPA 
are the: 
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• determination of the burden of pollutants (e.g. by using vehicle emission 
factors); 

• modeling of the dispersion of the pollutants around the source; 

• exposure assessment to evaluate the risk of the population exposed to the 
defined burdens; 

• evaluation of the impacts caused by the pollutants to the human health 
and to the environment; 

• monetary quantification of each impact (this step is usually based on the 
willingness to pay for reduced health risks). 

Costs of Up and Downstream Processes 

This costs category considers the external costs generated by indirect effects 
(not related to the transport journey itself) such as the production of energy, 
vehicles and transport infrastructure. These costs occur also in other markets, such 
as the energy market, in addition to the transport one so it is important to consider 
the appropriate level of internalization within these markets. The most relevant 
cost categories considered in [4] are: 

• energy production (well-to-tank emissions — WTT); 

• vehicle production, maintenance and disposal; 

• infrastructure construction, maintenance and disposal. 

The methodology adopted to calculate these costs is basically founded on the 
air pollution and climate change external costs estimation. The various studies 
treating this argument differ, among them, from the cost categories covered: for 
example some studies consider only the climate change costs of the up- and 
downstream processes whereas others also considers costs related to the air 
pollution costs. 



Noise Costs 

Exposure to noise emissions from traffic causes not only disturbs to people but it 
can affect their quality of life and health. The greater urbanization and the increase 
in traffic volumes are increasing the noise emissions. The two major negative 
impacts associated to this externality are: 

• costs of annoyance: due to social disturbances of persons exposed to 
traffic noise which result in social and economic costs such as discomfort 
and pain suffering; 

• health costs: noise level above 85 [dB(A)j can cause hearing damage 
while lower level (above 60 [dB(A)j may result in changing of heart beat 
frequency, increasing of blood pressure and hormonal changes. 
Furthermore, noise exposure can increase the risk of cardiovascular 
diseases and decrease the quality of sleep. 
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The most used methodology to estimate the marginal noise external costs is 

the Impact Pathway Approach. The key steps of the IPA [11] are: 

• level of noise emissions measured in terms of change in time, location, 
frequency, level and source of noise; 

• noise dispersion models used to estimate the changes in the exposure to 
noise according geographical locations, and measured in dB(A) and noise 
level indication; 

• Exposure -Response Functions (ERFs) showing a relationship between 
decibel levels and negative impacts of the noise; 

• economic valuation techniques of the negative impacts of the noise 
identified; 

• overall assessment to identify aggregated economic values taking into 
account all the negative impacts identified. 

Other External Costs 

The researches on the external costs, generally, focus only on the most 
important cost categories costs (such as air pollution costs, noise costs, climate 
costs or accident) neglecting other external costs categories. The reasons are 
mainly due to the complexity in the impact pattern and uncertainty in the valuation 
approaches. Methodologies for the calculation of these external costs are present 
only in few studies thus, presently, are not as sophisticated as for the most studied 
external costs categories [4]. The other external costs categories estimation 
considered are: costs for nature and landscape, cost to ensure water and soil 
quality, costs to ensure biodiversity losses, cost in urban areas (such as separation 
costs for pedestrian and costs of scarcity for non-motorized traffic). 



1.3 A Case Study from Automotive Industry Logistics 

Worldwide, 65,462,496 passenger cars have been assembled and 62,786,169 
registered (or sold) in 2013. Fig. 1.3 shows the worldwide passenger cars 
production from 2000 to 2013 [12] by macro-areas while Fig. 1.4 shows the 
worldwide passenger cars registrations (or sales) from 2005 to 2013 [13]. 
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Fig. 1.3 New passenger cars assembled worldwide from 2000 to 2013 [own Figure based on 

[12] data] 




Fig. 1.4 New passenger cars registered (or sold) worldwide from 2005 to 2013 [own Figure 

based on [13] data] 



A widespread analysis covering the Europe has been conducted to identify the 
flows of the new passenger cars from the assembly plants to each national 
European market. The main steps followed to perform the analysis have been: 

• Passenger cars definition. According to Eurostat [14] and OICA [15], 
a passenger car is a motor vehicle with at least four wheels used for 
the transport of maximum nine passengers including the driver. 
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• Identification of automotive assembly plants in Europe. More than 
100 passenger cars assembly plants have been identified in the Europe 
excluding Russia [16]. 

• Identification of new passenger car models assembled in each plant. A 
web-research performed on the Original Equipment Manufacturers 
official websites allowed identifying the car models assembled in each 
plant. 

• Collection of new passenger cars registration statistics in 2013 for the 
European Countries. These statistics have been found in public 
databases and, for some Countries, contacting privately national 
statistics associations. The statistics are split by passenger cars model. 

• Quantification of the new passenger cars flows from the assembly 
plants to the selected countries. The passenger cars distribution flows 
have been calculated by crossing the data collected in the previous 
steps. Each model assembled only in one plant (for example the Dacia 
Logan assembled only in the Romanian plant of Colibasi) provided 
quite certain information about the related distribution flows origin. In 
some cases, in 2013, a car model has been assembled in more than one 
plant (for example Audi A4 assembled in the German plants of 
Ingolstadt and Neckarsulm): it has been made the assumption that in 
each European country the 50 % of the Audi A4 flow came from 
Ingolstadt and 50 % from Neckarsulm. More accurate hypotheses 
have been made if available the total number of the cars assembled in 
the plants split by model. For example Opel Astra, in 2013, has been 
assembled in Gliwice (Poland), Riisselsheim (Germany) and Bochum 
(Germany) plants, respectively 100,886, 58,547 and 16,339 units. The 
flows in this case have been split in proportion to the number of cars 
assembled in these three plants. 



1.3.1 Inland Waterway Transport (IWT) 

The European Union is characterized by a network of inland waterways of more 
than 40,000 [km] [1], 29,172 [km] of which have been earmarked by 
Governments as waterways of international importance [17]. The most important 
European waterways are located in the South-East corridor, East-West corridor, 
Rhine corridor and North-South corridor (Fig. 1.5). As aforementioned, inland 
waterway transport accounted in the EU28, in 2012, only for 4 % of freight 
transport based on tonne-kilometres [t-km] [1]. 
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EAST-WEST < 




Fig. 1.5 : Overview of European inland Waterways [18] 

Analyses conducted for the “business as usual” outlook for 2040 [19] show that 
the modal share of the Inland Waterway Transport (IWT) in Europe will not 
increase significantly. The high capacity for bulk transport offered by IWT is 
presently mainly exploited by the agricultural, metal and petroleum industry and 
in Western Europe by the hinterland transport of maritime containers. 
Opportunities to raise the modal share of the IWT could be found in the following 
market not widely exploited by this transport modality such as: 

• paper industry; 

• automotive industry; 

• high and heavy transport (e.g.: construction equipment, generators, 
turbines, wind turbines, etc.); 

• waste industry. 

One sector of the automotive industry with a huge potentiality for the IWT is the 
transport of the new passenger cars from the assembly plants to final customers. 
This market is, currently, not widely exploited in Europe where, however, since 
1997, one of the largest European automobile logistics service providers, 
transports new passenger cars on the Danube waterway. New passenger cars of 
Ford and Renault are transported on the Danube from Kelheim to Budapest and on 
the way back Suzuki cars are transported from Budapest to Kelheim. Every year, 
Suzuki transports about 18,000 vehicles on the Danube according to a regular 
schedule with two departures per week by two self-propelled vessels. The 
“Kelheim” and the “Heilbronn” (see Fig. 1.6) can load up 200-260 cars on three 
decks. The loading and unloading of the vehicles are carried out through a ramp 
installed at the bow of the ship, which can be lowered on the Ro-Ro ramp of the 
port [20]. 
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Fig. 1.6: Heilbronn vessel [21] 



As aforementioned, more than 100 passenger car assembly plants were present 
in the EU28 in 2013 [16]: among these, 37 were present in ten Countries of the 
Rhine-Maine-Danube area considered in [22] (Fig. 1.7). 




Fig. 1.7 : Potential Countries for the distribution of the new passenger cars in the Rhine-Main- 

Danube area [22] 



Table 1.1 shows the quantification of the potential passenger car flows among 
these Countries. 
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1.3.2 Discussion 

A numerical example tried to quantify the potential external costs savings due to 
the 1WT to transfer new passenger cars from four selected assembly plants located 
in Romania (2), Hungary (1) and Slovakia (1) to Germany, Netherlands and 
Belgium [23]. It has been assumed to unload: (1) in the Port of Rotterdam all the 
passenger cars directed to Belgium and Netherlands; (2) in the Port of Mainz the 
40 % of the new cars directed to Germany. The Marco Polo Calculator has been 
used for calculating the external costs of the transport. The European Union’s 
Marco Polo Programme [24] aims at shifting or avoiding freight transport off the 
road to other more environmentally friendly transport modes. This Programme, 
running by yearly calls for proposals, selects for financial support the proposals 
received depending on the level of the environmental and social benefits expected 
by them. The Marco Polo Calculator is a tool developed for the applicants of the 
Marco Polo Programme allowing comparing the monetized external cost impact. 
This calculator provides external costs estimates for road, rail, inland waterway 
and short sea shipping transport modes. In Table 1.2 the external costs coefficient 
used for the case study are shown [25]. 



Table 1.2 External costs of the freight transport [€cent/t km] as in Marco Polo Calculator 



Categories 


Road 

(motorways) 


Electric 

Train 


IWT 

Capacity: 401-650 [t] 
Fuel: LNG 


Air Pollution 


0.858 


0.100 


0.180 


Climate Change 


0.392 


0.146 


0.120 


Noise 


0.193 


0.149 


- 


Accidents 


0.064 


0.033 


- 


Congestion 


0.343 


0.020 


- 


Total 


1.85 


0.448 


0.3 



Vessels propelled by Liquefied Natural Gas (LNG) have been chosen to perform 
the case study among the many fuel technologies available for the inland 
waterway vessels. Table 1.3 shows an estimate of the reduction of CCL and air 
pollutants (PM, SO 2 and NO x ) emissions comparing low sulphur fuel oil and four 
alternative fuel technologies [25]. 

Table 1.4 shows the flows of the new passenger cars, in 2013, from four selected 
assembly plants to different national markets (Belgium, Netherlands and 
Germany). 

Table 1.5 shows, instead, the external costs reduction of the multimodal transport 
(with road transport or rail transport as main haulage) compared to the direct road 
transport. In order to facilitate the comparison, post haulage has been neglected. 
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Road and rail distances, for a given pre-haulage route, have been assumed to be 
equal and calculated using Google Maps [26] instead Ecotransit calculation tool 
[27] provided the river distances for the main-haulage route. 

Despite the river distance always more than 20 % longer respect to the direct road 
distance for all the selected routes, the multimodal transport with IWT as main 
haulage and rail transport as pre -haulage is the best solutions in terms of total 
external costs resulting in a cut of around the 80 % compared to direct road 
transport. 



Table 1.3 CO2 and air pollutants (PM, SO2 and NO x ) reduction for inland vessels considering 

different fuel technologies 



Fuel technology NO x PM SO 2 CO 2 Fuel Consumption 



Fuel oil (low sulphur) 
Diesel Particulate filter 
(DPF) 


— 


-68 % 


Base option 

+2% 


Selective Catalytic 
Reduction (SCR) 


-85 % 


- 


- - 


DPF + SCR 


-85 % 


-68 % 


+2% 


Liquified Natural Gas 
(LNG) 


-75 % 


-97 % 


- -10 % 




Table 1.4 Number of new passenger cars transported (year 2013) 
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Table 1.5 External costs reduction of multimodal transport 
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Abstract 

The Economic Order Quantity (EOQ) model, proposed by Harris in 1913, is 
one of the most studied models for the inventory management. The model aims at 
identifying the optimal lot size minimizing the total inventory costs, typically only 
holding and ordering costs. Many researchers extended this model trying to adapt 
it to real-life situation by providing new mathematical models. The increasing 
attention paid to sustainable manufacturing led, in the last years, to include the 
external costs of logistics in the EOQ model. Starting from a literature review on 
the inventory management models, this chapter defines the new Sustainable Order 
Quantity (SOQ) model. In the model, the loss factor parameter quantifies the loss 
in energy per unitary load transported and unitary distance covered, univocally 
identifying the various transport means. The optimal order quantity is derived 
minimizing a logistic cost function that considers both economic and social- 
environmental costs. The model allows determining at the same time the reorder 
level, the safety stock as well as the optimal transport means. 

Keywords : Inventory management. Economic Order Quantity (EOQ), 
Sustainable Order Quantity (SOQ), Logistic cost function. Loss factor of 
transport. Environmental cost, Transport means selection 

2.1 Notations 

Notations adopted in this chapter are in Table 2. 1 . 



Table 2.1 Notations adopted 



Variable 


Name 


Unit 


/ 


Loss factor 


- 


G 


Product annual requirement 


[unit/year] 




Logistics costs 


[€/year] 


ftp 


Purchase costs 


[€/year] 




Ordering costs 


[€/year] 


d> H 


Holding costs 


[€/year] 




Transport costs 


[€/year] 


@S 


Shortage costs 


[€/year] 


&EX 


External costs 


[€/year] 


Q 


Order quantity 


[unit] 


r 


Reorder level 


[unit] 


H 


Production hours per year 


[h/year] 
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m 


Mass of one product 


[kg/unit] 


A 


Product demand of the z'-th period 


[unit/period] 


E(Di) 


Expected value of the product demand in the z-th 
period 


[unit/period] 


&Di 


Standard deviation of the demand in the z-th period 


[unit/period] 


LT 


Supply lead time 


[h] 


E(LT) 


Expected value of the supply lead time 


[h] 


&LT 


Standard deviation of the supply lead time 


[h] 


L 


Transport distance 


[km] 


CT 


Consumption time 


[h] 


c H 


Unit annual holding cost 


[€/unit-year] 


Co 


Order cost 


[€/order] 


SOQ 


Sustainable order quantity 


[unit] 


SS 


Safety stock level 


[unit] 


D* 


Maximum demand value of the z-th period not 
causing a stock-out event in that period 


[unit/period] 


Dtot 


Demand on LT periods 


[unit/period] 


D tot 


Maximum value of the demand on LT periods 
not causing a stock-out event 


[unit/period] 


I 


Inventory level 


[unit] 


LT * 


Maximum value of the supply lead time not 
causing a stock-out event at a given service level 


[h] 


N s 


Number of shortages in one ordering cycle 


[unit] 


cs 


Unitary shortage cost 


[€/unit-order] 


C T 


Unitary transport cost 


[€/t] 


(a,b,c) 


Regression parameters of transport costs functions 


[€/kg] 


(ki,k 2 ,k 3 ) 


Regression parameters of average speed functions 


[km/h] 


Vact 


Transport means actual speed 


[km/h] 


V 


Transport means average speed 


[km/h] 


T t 


Transport time 


[h] 


T l 


Time required for the material handling, order 
management and quality control 


[h] 


k 


Ratio between T t and T L 


- 


Er 


Energy consumption per order required by a given 
transport means 


[MJ/order] 


Si 


Monetary cost per unit mass emission of the z-th 
pollutant 


[€/kg] 


£ex 


Unit cost of externalities 


[€/t-km] 


E s 


Energy consumption per functional unit 
transported 


[MJ/tkm] 


(aJXYi) 


Regression parameters of transport emissions 
functions 




(S,M) 


Regression parameters of unit external costs 
functions 
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2.2 Overview of the State of the Art 

The Economic Order Quantity (EOQ) model [1] is one of the most investigated 
inventory management model aiming at identifying for a single inventory item its 
optimal lot size under the following hypothesis: 

(i) product demand constant over time; 

(ii) negligible supply lead time; 

(iii) unit production cost independent from the production quantity; 

(iv) zero-defective products; 

(v) fixed transport costs (implicitly considered in the ordering costs). 

Under these assumptions, the optimal order quantity is obtained as a trade-off 
between holding and order costs (Eq. 2.1). 

EOQ = \ 2 ' G ' C ° (2.1) 

V C H 



Starting from the work of Harris, many models have been developed including 
other cost figures in the optimization function. 

A first class of models explicitly considers transport costs in the EOQ model. 
In logistic systems these costs depend on the transport means adopted, as well as 
on the shipment size [2], and could weight upwards of 50 % of the total logistic 
costs. The transport costs are explicitly considered in [3]; these costs are also 
considered in [4] by adopting freight rate functions available in the scientific 
literature. In [2] an optimal solution procedure for solving the EOQ models is 
provided in case of transport costs are explicitly considered and shaped as all-unit- 
discount costs. In [5] optimum lot-sizing algorithms in case of quantity and freight 
discounts are proposed: all-units and incremental discounts are considered. 

In a further class of models restrictive hypotheses of [1] are released. In [6] 
product quality problem is investigated, and the replacing of a random quantity of 
defective items is considered under an EOQ model. The effects on the inventory 
management of imperfect quality items are considered in [7]; in the model, the 
percentage of the imperfect items is characterized by a known probability density 
function assuming also that they can be used in another production/inventory 
context, generating less revenue than good quality ones. The model proposed in 
[8] extends this model taking into account the probability of failures in inspection 
activities. In [9, 10] an EOQ model under the hypothesis of exponentially 
decaying inventory and constant product demand is defined. Authors in [11] 
propose a model allowing evaluating the EOQ value in case of deteriorating items 
and permissible supply delay. 
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The EOQ model proposed in [12] quantifies the increase in the optimal order 
quantity values due to the inflation effects on the prices; the effects of inflation 
uncertainty on inventory decisions are also investigated in [13]. Models including 
simultaneously inflation and deteriorating items effects are proposed in [14] and 

[15]. 

In [16, 17] the assumption of a deterministic demand is released, and the EOQ 
model is extended considering, respectively, the case of a price dependent demand 
and of a demand varying stochastically. Many contributions are available in which 
traditional EOQ model is applied under the hypothesis of stochastic variability of 
product demand; in case of uncertain product demand the safety stock sizing 
problem is integrated in the EOQ model. In [18] a model allowing obtaining 
optimal service level and safety stock level as a function of the ratio Q/o x is 
proposed, where o x is the standard deviation of the supply lead time. The 
evaluation of the optimal safety stock levels, through the EOQ model, of 
components assembled to obtain a finished product is proposed in [19]. The 
effects of partial backordering on EOQ solution in case of a variable product 
demand is evaluated in [20] by means of a deterministic inventory model: during 
stock out periods a fraction of the demand is backordered and the remaining 
fraction generates shortage costs. In [21] a solution procedure to compute EOQ in 
case of backordering is provided analyzing two different optimization problems, 
providing the optimal value of the maximum inventory level and the optimal value 
of the backorder level, respectively. The analysis of the potential benefits of 
Vendor-Managed Inventory (VMI) implementation in EOQ model is proposed in 
[ 22 ]. 

Lead time variability is investigated in the model proposed in [23]: lead time is 
assumed as a decision variable obtaining its optimal value by means of 
minimizing crashing costs, defined as extra costs to be charged in order to reduce 
lead time. In [24] optimal lead time value, as well as order quantity and safety 
stock values are obtained in case of crashing lead time costs and price discounts of 
backorders. In [25] an algorithm aiming at solving the single vendor single buyer 
problem in case of stochastic variability of lead time demand and a lead time 
linearly varying is proposed. The EOQ model is applied in [26] considering 
stochastic variability of lead time and deterministic demand rate. In [27] optimal 
order quantity and reorder level values are obtained in case of random lead times 
by assuming that it is possible to obtain expediting orders (orders with a shorter- 
than-average lead time at an extra cost). Methods to reduce supply lead time 
variability are analyzed in [28]: order splitting is identified as the optimal solution 
in case of a lot size-dependent supply lead time. The effect of reducing lead time 
variability on safety stock level, in case of a gamma distributed lead time, is 
investigated in [29]. Authors in [30] show how the reduction of lead time 
variability is more effective than the reduction of its expected value in case of a 
deterministic product demand. 
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Sustainable inventory management is a quite new research field. The 
increasing attention paid to sustainable manufacturing led to the development of 
logistic models aiming at jointly minimizing internal and external costs of 
logistics. Many models are available in scientific literature mainly focusing on the 
reduction of the carbon footprint of logistic activities. Less frequently the other 
categories of the external costs are considered. Authors in [31] studied the 
environmental impact in the inventory problem considering the carbon emissions 
due to the energy and fuel consumption in the transport and storage of 
perishable goods. In [32] authors consider an environmental performance based 
green cost as a linear function of the production/order quantity in the economic 
production quantity (EPQ) model and in the economic order quantity (EOQ) 
model; results of a sensitivity analysis showed that for the two models the optimal 
quantity value is smaller than in the standard case when green costs are accounted. 
In [33] social and environmental costs have been included into the classical EOQ 
model. Authors included social aspects in term of working hours and considered 
five environmental approaches to account the carbon footprint: direct accounting 
(carbon footprint considered as an additional economic cost), carbon tax (applied 
by regulatory agencies), direct cap (imposed by regulatory agencies or 
alternatively by public awareness about green products), cap&trade (rewards for 
company emitting less than an allowed cap and penalization in the other case) and 
carbon offsets (investments made to reduce emissions such as efficient and 
renewable energy resources). Similar considerations are in [34], where authors 
show how the overall emissions across the supply chain can be significantly 
reduced in case of collaborative firms within the same supply chain. In [35] 
authors examined the impact of the emission trading mechanism in the inventory 
management. They considered greenhouse gas emissions due to transport and 
warehousing operations. In [36] authors presented an overview about the 
inventory management aspects causing environmental damages that are not 
analyzed in depth by the traditional inventory analysis such as packaging, waste 
and location of stores and proposed a model that includes into the classical EOQ 
model environmental impacts related to vehicles journeys and waste disposal. In 
[37] a two-level supply chain model is proposed to determine the optimal 
shipment size and the optimal number of shipments taking into account the CO 2 
emission costs due to transport operations. The environmental costs are split into 
two terms: fixed (depending on the vehicle type, vehicle age and average speed) 
and variable costs (depending on the actual weight of the shipment). In [38] the 
classical EOQ model is extended in order to solve a multi-objective problem 
where the goal to be achieved is not only the minimization of the logistics costs, 
but also the reaching of environmental and social goals. In [39] authors provided 
actual external costs evaluations based on a wide analysis of international 
researches carried out on this topic; 13 categories of external costs generated by 
the transport means are considered and a sensitivity analysis shows how the 
distribution network configuration could be influenced by a variation of these 
costs. In [40] authors proposed a model aiming at extending the external costs 
analysis of the inventory management; results of a life cycle assessment of the 
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order quantity supplied are presented; environmental and social costs related to the 
transport, handling, storage, and waste disposal are considered. 

In [41] authors proposed a Sustainable Order Quantity allowing evaluating the 
sustainable order quantity jointly minimizing the logistic and the environmental 
costs, as well as to identify the proper transport means under deterministic demand 
and lead time; the SOQ model jointly solves the two problems by the adoption of 
the “loss factor”, a parameter which measures the loss in energy occurring during 
the shipping of materials by a given means of transport [42]. The global warming, 
the acidification and the tropospheric impacts are considered in the external costs 
function. In [43], authors adopted model in [41] in order to investigate the effects 
of different internalization strategies of external costs of transport on the optimal 
solution of the logistics problem. In [44] the model has been extended by 
considering a stochastic variability of the product demand. A further extension of 
the model in [41] is proposed in [45] by adding the environmental costs due to the 
spare parts repairing and purchasing. In [46] uncertain demand of repairable spare 
parts is considered. In [47] a stochastically variable lead time is considered, and 
the external costs due to the freight transport means as defined in European 
Commission commissioned studies on this issue are evaluate. 



2.3 The Loss Factor of Transport 

It is well known that if we want to move a load, we need to spend a certain 
amount of energy. This energy is partially transferred to the load, and partially lost 
due to the friction: the first amount could be theoretically regained but the latter 
have to be considered a pure loss, since there is no way to recover it. 

Considering the simple case of a load W pulled over a horizontal floor (see Fig. 
2.1a), it is possible to calculate the driving force F needed to move the load as: 

F = f-W, ( 2 . 2 ) 

with / the coefficient of friction between the lower surface of the load and the 
floor. In the more general case, the total amount of energy required for transport 
the load W along a route of length L is: 

E = A E p + A E k + F L, (2.3) 



where 

AE P and AE K are the differences between the initial and the final potential and 
kinetic energy of the load, respectively. In case of these differences are negligible, 
the amount of energy required is equal to F L. Being F = f W, it results 

E = f -W ■ L. (2.4) 

When the load W is shifted for the distance L using two wheels of negligible 
weight (see Fig. 2.1b) a smaller amount of energy {E') is required, due to the 
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different value assumed by the frictional coefficient: Coulomb frictional factor in 
the first case (/), a rolling frictional factor in the second one (/'): 

E' - f'-W-L. (2.5) 

It is not possible to easily apply Eq. (2.3) in real transport system in order to 
compute the exact amount of energy required for the transport, but it is possible to 
measure the amount of energy E required to transport a given load W along a route 
of given length L. 

The amount of energy required E depends on the transport means adopted. The 
ratio E/(W’L), has been defined in [42] as the loss factor of transport (/). The loss 
factor firstly appeared in scientific literature in [42] as a technical energy-related 
performance measure conceived for taxonomy purposes. By adopting available 
data on freight transport systems, the author [42] provides a / taxonomy of both 
continuous and discontinuous transport means (see Fig. 2.2). According to the 
Jonkers taxonomy, a given means of transport is characterized by an / value, 
which is an average measure for a given type of transport means. 
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Fig. 2.1 Different ways of transporting a load 
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Fig. 2.2 Loss factors of different transport means (source: [42]) 



Loss factor values are strictly related to the energy conversion technology of 
transport means. In the last years, these conversion technologies significantly 
evolved due to technical improvements (e.g. new fuel injection technologies), as 
well as to policy regulations (e.g. fuel composition, emission limits). In Table 2.2 
loss factor values evaluated from 2009 data of [48] database are shown. In Table 
2.3 data from [48] (year 2012) and [49] (year 2012) are listed. 

Table 2.2 Loss factor (f) value for different means of transport (data 2009) 



Transportat means 


f 

m 


Ship — Oversee 


0.010 


Ship — Coast 


0.033 


Rail — Diesel 


0.033 


Rail — Electric 


0.037 


Ship — Inland 


0.048 


Truck" 32 - 40 [t] 


0.118 


Truck" 28 - 32 [t] 


0.134 


Truck 20 - 28 [t] 


0.162 


Truck" < 28 [t] 


0.170 


Truck 14 - 20 [t] 


0.214 


Truck 7.5 - 14 [t] 


0.242 


Truck < 7.5 [t] 


0.475 


Small truck <3.5 [t] 


0.718 


Airplane — international 0.989 


Airplane — inland 


1.774 



a Operated with trailer 
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Table 2.3 Loss factor (f) value for different means of transport (data 2012) 



Transport means 


/ 

J/J 


Ship — Oversee 


0.010 


Ship — Tank 


0.021 


Ship — Coast 


0.033 


Rail — Diesel 


0.033 


Rail — Electric 


0.037 


Ship — Inland 


0.048 


Truck 34 - 40 [t] 


0.063 


Truck 28 - 34 [t] 


0.069 


Truck 26 - 28 [t] 


0.119 


Truck 20 - 26 [t] 


0.136 


Truck 14-20 [t] 


0.127 


Truck 12- 14 [t] 


0.139 


Truck 7.5 - 12 [t] 


0.149 


Truck 3.5 -7.5 [t] 


0.165 


Truck 3.5 -7.5 [t] a 


0.156 


Truck <3.5 [t] b 


0.335 


Truck <3.5 [t] c 


0.554 


Airplane 


0.984 



d Truck operated without trailer 



b Light good vehicles (LGV) 
“Vans 



2.4 A Sustainable Order Quantity (SOQ) Model 

The inventory management model proposed in this work extends the definition of 
the Economic Order Quantity (EOQ) and defines a Sustainable Order Quantity 
(SOQ) as the order lot size jointly minimizing a logistics cost function @ L 
considering both logistics and external costs. By identifying with a proper loss 
factor / each transport means, it follows 

min <P L [€/year], (2.6) 

The annual logistics cost function <P L is defined as the sum of the purchase, 
ordering, holding, transport, shortage, repair, and external costs. All costs are 
evaluated on a yearly basis. 

+ 0 0 + + O r + [€/year] (2.7) 



In the following, each cost figure of (2.7) is detailed. 




30 2 Sustainable Inventory Management 



2.4.1 Purchase and Ordering Costs 

With G the annual product requirements, in case no quantity discount are 
considered, the purchase costs of the items can be computed as the product of the 
annual amount ( G ) and the unit purchase cost. In absence of quantity discount, 
however, purchase costs do not affect the solution of the logistics problem. 

With G/Q the average number of ordering cycles in 1 year, the ordering cost is 
evaluated as: 



Q 

®o= c o-— [€/year] 



(2.8) 



2.4.2 Transport Costs 

Transport cost is the sum of the material handling and of the shipping costs. Many 
models are available in scientific literature allowing shaping the transport cost as a 
function of the transport distance for an assigned class of transport means [4, 50, 
51]. For a given transport means and transport distance, the cost of transport 
depends on both the shipment size [2] and the speed of the transport means [52, 
53]. In the model proposed, each transport means is identify by a proper loss 
factor value, which represents an average value of its specific energy 
consumption. As a consequence, for each means of transport an average 
characteristics speed value is considered. A single-modal transport is considered 
in the model. Multimodal option can be analyzed as an extension of the model 
proposed by considering the overall route length as the sequence of segments each 
operated by a different transport means. Finally, partially loaded transport means 
configurations are not considered. 

Under these assumptions, the unit transport cost ( c T [€/t]) can be considered a 
function of the transport means adopted and of the transport distance: 

c T = c T (f ,L) [€/t], (2.9) 

and the yearly transport cost can be obtained as: 

O T =G-c T (f,L) [€/year]. (2.10) 



2.4.3 Holding Costs 

The holding costs are evaluated by means of different functions in the three cases 
considered: 

(a) deterministic demand and lead time; 

(b) stochastic demand and constant lead time; 

(c) deterministic demand and stochastic lead time. 
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(a) Deterministic demand and lead time 

In case of a supply lead time and a product demand constant over time (see Fig. 
2.3), the holding costs can be evaluated as: 

®h=Y c h'Q [€/year] (2.11) 




Fig. 2.3 Inventory level (I) over time in case of constant demand and lead time 



(b) Stochastic demand and deterministic lead time 

Stochastic variability of the product demand leads to stock-out events at a given 
service level ( SL ) (i.e. probability of no stock out). The adoption of safety stocks 
allows facing with stock-out events. In case of a safety stock is adopted, the 
holding costs are: 



^ H ~ C H 



— Q + SS 
2 



[€/year] 



( 2 . 12 ) 



In Fig. 2.4, order quantity Q, safety stock level SS and reorder level r are depicted. 
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Fig. 2.4 Inventory level (!) over time in case of stochastic demand 



The required safety stock level can be evaluated starting from the service level and 
the demand variability provided by the production system. With D, the product 
demand of the z-th period of the lead time (LT), under the following assumptions: 

- Dj are independent stochastic variables; 

- product demand is characterized by the same expected value and by the 
same standard deviation in each z-th period of the lead time: 

E(D,) = E(DJ = ... = E(D lt ) = E(D) = G/H, 



2 2 

(T n (7 n . . . ■ 



-Or 



= 07 



The safety stock level, SS, consistent with an assigned service level, SL , can be 
computed as: 



SS = LT ■ 



D ' -£(£>)] 



[unit]. 



(2.13) 



The service level (SL) is: 

LTD ' 

SL = prob | D tot <LT-D* j = J pdf(D TOT )dD TOT , (2.14) 

—oo 

with D* the maximum demand value of the z-th period not causing a stock-out in 
that period, D T ot the demand on LT periods, and LTD* = D T or the maximum 
value of the demand on LT periods not generating a stock-out (see Fig. 2.4). 

(c) Deterministic demand and stochastic lead time 

In case of a stochastic variability of supply lead time, the expected inventory level, 
E(1 ), in one ordering cycle have to be computed. The inventory level ( L) is affected 
by both the order quantity size ( Q ) and the safety stock level ( SS) as well as by the 
supply lead time variability. 



2.4 A Sustainable Order Quantity (SOQ) Model 33 



The safety stock level value ( SS) consistent with an assigned service level (SL) can 
be evaluated as: 



SS = D [LT* -E(LT)] [unit]. 
The service level (SL) is: 



(2.15) 



LT 



SL = prob(D TOT <LT* D) = J pdf(LT)dLT = F(LT*), (2.16) 



with LT * the maximum value of the supply lead time not causing a stock out event 
at a given service level (SL) (see Fig. 2.5), and D TO t the lead time demand. 
Introducing the standardized variable (z) of the stochastic variable LT, it is 
obtained: 

SS = D ■ z* ■ < j lt [unit], (2.17) 



SL = F(z*), 



with 



LT* -E(LT) 



(2.18) 



(2.19) 




Fig. 2.5 Inventory level (I) in case of stochastic lead time ( LT) and LT= E(LT) 



Three different situations may occur due to the stochastic variability of LT: 

(A) LT <E(LT) (Fig. 2.6), 

(B) E(LT)<LT<LT*( Fig. 2.7), 

(C) LT>LT * (Fig. 2.8). 
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Fig. 2.6 Inventory level (L) in case of stochastic lead time (LT) and L T < E(LT) 




LT * 



Fig. 2.7 Inventory level (I) in case of stochastic lead time ( LT) and E(LT) < LT < LT* 




Fig. 2.8 Inventory level (1) in case of stochastic lead time (LT) and LT> LT* 
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The expected inventory level during one ordering cycle, E(I), can be evaluated as: 



In the reference case of LT = E{LT) (see Fig. 2.5, and the thin line in Fig. 2.6), 
the average inventory level and the consumption time are respectively equal to: 
E(I) = 1/2 Q + SS, and CT = Q/D. In cases B and C the average inventory level is 
affected by the lead time variability in the next ordering cycle; in case C, the 
length of the ordering cycle is greater than Q/D. In Table 2.4 the expected 
inventory level, the consumption time, and the corresponding occurrence 
probability values in the three cases considered are illustrated. 

In case of high service level values (SL > 0.90) the following considerations arise: 

• the term [LT-LT ] assumes small values; 

• the occurrence probability of case C is very small; 

• the consumption time in case C can be approximated to Q/D; 

• the expected inventory level in case C is equal to ‘/i-Q; 

• the expected inventory level E(l) in one ordering cycle is obtained as: 



fE(LT) 

E(I) = E (I) A ■ pdf (LT)dLT 




[unit/order]. 



(2.20) 




[unit/order], (2.21) 



and the corresponding holding cost is: 




[ pdf (z*) + z * -F(z*)] 
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Table 2.4 Expected inventory level and ordering cycle length in the three cases considered 



Expected inventory level E(I) 


Case A 
LT < E(LT) 


E(I a ) = \q + D ■ [E(LT) - LT] + SS 


Case B 

E(LT)<LT<LT* 


E(I b ) = ^Q-D[LT-E(LT)] + SS 


Case C 
LT> LT* 


e{i c ) = \q+\ss 


1 [LT*-E(LT)\ 

2^ CT c 




Consumption time CT 


Case probability 


Case A 
LT < E(LT) 


Q 

D 


\pdf(z)dz 

—CO 


Case B 

E(LT)<LT<LT* 


Q 

D 


\pdf(z)dz 

0 


Case C 
LT> Uf 


—+ILT — LT*] 
D 


+00 

J pdf ( z)dz = l-SL 

z* 



2.4.4 Shortage Costs 

In case of stock out events occur, shortage costs have to be considered. They can 
be evaluated for both the two cases: 

(a) Stochastic demand and deterministic lead time 

(b) Deterministic demand and stochastic lead time 



(a) Stochastic demand and deterministic lead time 

The number of shortage units per ordering cycle, N s , can be calculated as: 






tCXj 

j* pdf (D tot ) (d-^TOT — ^4 tot ) dDjQj [unit] . 



Under the assumptions in Sect. 2.4.3, since 

D tot = D x + D 2 + ... + D lt [unit/h] 

it follows: 

E(D tot ) - LT ■ E(D), 



(2.23) 

(2.24) 

(2.25) 



2 t rji 2 

0"d J^l ' (Tjj. 

u TOT u 



(2.26) 
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For any pdf{D ,), as stated by the central limit theorem, the pdJ{D TO T ) tends 
towards a normal distribution: 



pdf(D TOT ) — 



J~LT - a n • V2tt 



■exp 



(D tot -LT -E (D)) 2 



2-LT-al 



. (2.27) 



The yearly shortage cost, 0 S , can be evaluated as: 

G 



^ s Qs ’ ’ f P4f (.D tot ) -{D T ot D T oT)dD T oT [€/ycsr]. (2.28) 

Q D * 



(b) Deterministic demand and stochastic lead time 

Stock-out events occur only in case C ( LT > LT *) discussed in Sect. 2.4.3, as 
shown in Fig. 2.8. The average number of shortage units per ordering cycle can be 
evaluated as: 



-l-rxj 

N S = D ■ J (LT -LT' )■ pdf(LT) ■ dLT [unit/order], (2.29) 



and the corresponding annual shortage cost as: 

d ) s =c s ~-D-a LT -L{z*) [€/year], 

with L(z*) the standardized normal loss function defined as: 



(2.30) 



L(z*)= ^ (z - z*) ■ pdf (z) dz = /?^(z*)-z*-[l-F(z*)]. (2.31) 



2.4.5 External Costs 

In the logistic cost function of the SOQ model (Eq. 2.7) costs of externalities are 
considered. As discussed in [40], different inventory management activities are 
responsible of environmental and social costs. The main contribution in terns of 
externalities has to be addressed to transport activities. Transports are responsible 
of environmental impacts, as well as of social damages. 

As far as concern environmental impacts of transport, the related external costs 
can be computed by adopting two different approaches. 

In [41, 44-46] the environmental costs (@ E [€/year]) for a given freight 
transport means are modelled as a function of the loss factor as: 
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0 



E 



Qi E s w 



[€/year], 



(2.32) 



where: 

E s [MJ/tkm] is the energy consumption per functional unit transported; 

E R (f) =E s -Q-m L [MJ/order] is the energy consumption per order required by 
a given transport means, being m the mass of the single unit transported; 

Cj(f) [kgi,eq/t-km] are the environmental impacts (i = Global Warming, GW, 
Acidification Potential, AP, and Tropospheric Ozone Precursor Potential, TOPP 1 ) 
of the emissions of the transport means adopted, expressed as mass of the 
equivalent /-th pollutant (/ = (C0 2 ) eq ., (S0 2 ) eq ., NMVOC eq .) emitted per functional 
unit [tkrn]; 

e, [€/kg] is the unit cost per impact category. 

Environmental impacts are evaluated as: 

e,CO = 2XA (/) [kgi.eq/t • km]- (2.33) 

k 

where bk(f) [kgi/t-km] is the emission factors of the k - th pollutant of a given means 
of transport, and c,k [kgj eq / kg k | is the relative contribution (characterization 
factor) of the k - th burden to the i - th impact category. 

An alternative approach is the one adopted in [43, 47], where estimates of 
external costs of freight transport available in scientific literature [54, 55], and in 
EU official guidelines [56, 57] are assumed. 



1 Tropospheric Ozone Precursor Potential; chemical compounds, such as CO, CH 4 , NO x , and Non- 
Methane Volatile Organic Compound (NMVOC), reacting in the presence of solar radiation and 
producing ozone in the troposphere; they are expressed in NMVOC equivalents, weighed factors are in 
[58]. 
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3 SOQ Model Formulations 



Abstract 

This Chapter, starting from the definition of the logistic cost function defined in 
Chap. 2, presents the Sustainable Order Quantity (SOQ) model formulation 
considering product demand and supply lead time firstly as deterministically 
known and then characterized by stochastic variability. In the stochastic problem, 
solution procedures are suggested for solving the SOQ problem. Five cases 
(numerical experiments and industrial full-scale) are presented to test the model. 
Sensitivity analyses are carried out varying costs, distances, and uncertainty of 
demand and lead time, thus showing the model capabilities and limitations. 
Finally, results are discussed and compared with classical EOQ solutions. 

Keywords : Sustainable Order Quantity (SOQ), Stochastic variability of 
product demand, Stochastic variability of supply lead time, Sensitivity analysis, 
Automotive supply chain. Spare parts, Repair policy 

3.1 Deterministic Demand and Lead Time 

The following assumptions are assumed for the deterministic SOQ model 
formulation: 



1 . the annual demand is known; 

2. the item consumption rate is constant over the year and deterministically 
known (Q = G/H)\ 

3. transport is carried out by a single means of transport; 

4. the transport cost per order and the transport lead time depend on the transport 
means adopted for a given route; 

5. the supply lead time is calculated as the sum of transport time (TV) and the 
time required for the material handling, order management and quality control 
(7i); T l is assumed to be a fraction ( k < 1) of 7>; the transport speed is 
calculated from the free flow speed ( v ac , ) of different transport means as 
v = v ac ,/(l + k)\ the value of the supply lead time is: LT = L/v. 



In case of a deterministic demand and lead time, the total annual logistic cost 
function can be formulated as: 



_ Q G G 

^ L — C h ' 2 + Q ' C 0 + Q ‘ C T (f) + ^ EX 



[€/year]. 



(3.1) 



The optimal values of loss factor (fopp ) and of the sustainable order quantity 
(SOQ) jointly minimizing the logistic and external costs can be obtained by 
solving problem (3.2): 
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min<P,. (3.2) 

SOQ.f 

In case of deterministic demand and supply lead time, stock out events do not 
occur. For this reason, in (3.1) shortage costs are not considered. Under the 
assumptions made, order quantity ( Q ) can be obtained as: 



Q = 



GL 
pH v 



(3.3) 



Parameter p is the ratio between the reorder level (r) and the order quantity ( Q); 
it ranges in ]0; 1] interval. For a given value of the average transport speed, the 
corresponding minimum value of Q can be evaluated assuming p = 1 (r = Q). 
This value represents the minimum lot size avoiding the stock out when a means 
of transport (and consequently a given average transport speed) has been selected 
(see Fig. 3.1). 




Fig. 3.1 Lot size (0 vs. transport speed (v) 



By substituting Eq. (3.3) into Eq. (3.1), the annual logistic cost function is: 

[€/year], (3.4) 



®l=G- 



f r r \ 1 ^ 

m-c T (f ,L)+ — - 



2 p-H-v 



Starting from Eq. (3.4), a logistics, no-dimensional cost factor, F,, is defined by 
dividing the total annual costs by the constant value G-C/J2; this value represents 
the annual inventory cost occurring in a limited situation where the order quantity 
is equal to the annual requirement (G): 
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(3.5) 



The general logistic optimization problem is finally defined as: 



min F, . 

{ fp J 



(3.6) 



The reorder level is, usually, assumed by the logistics management based on 
the lead time required from the order emission to the inventory replenishment. In 
such circumstances, /is the only optimization variable in the problem (3.6) since p 
values are assigned. 



3.1.1 Environmental Costs 

As discussed in Sect. 2, there are many categories of externalities due to the 
transport. In this Section, the emissions of transport and its related costs are 
considered in the logistics cost function. In case environmental costs of transport 
are considered in Eq. (3.1), the logistic cost function is obtained by substituting 
Eq. (2.32) and Eg = E s Q m L in Eq. (3.5). 



Starting from data in [1, 2] on specific energy consumption (E s ), the 
corresponding loss factors values (/) and emission factors (e : ) have been evaluated. 
Values obtained are in Table 3.1, where estimates of the transport speed adopted 
are also listed. Data adopted here are expressed in [62009]- 




(3.7) 
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A regression analysis carried out on data in Table 3.2 led to shape the 
dependency of the average speed of transport (v), obtained with a k value of 0.5, 
and of impact factors ( e,) on loss factor (/) by means of quadratic functions: 



v = K-f 2 +k 2 -f + k 3 , 


(3.8) 


e i =r i -f 2 +a i -f +A- 


(3.9) 



Regression parameters as well as R 2 values are in Table 3.2. 

Table 3.2 Results of the regression analysis 



Impact category 


Unit 


Regression parameters 

cci Pi r 


R 2 


Global warming 


COjeq. 


0.8115 


0.0097 0 


0.9812 


Acidification 

Potential 


S02eq. 


0.0027 


0.0005 0 


0.9182 


Tropospheric 
Ozone Precursor 
Potential 


NMVOC er , 


0.0129 


0.0000 -0.0123 


0.9747 






ki 


k-2 kj 




Transport speed 


rkm/hl 


731.9 


-370.6 65.4 


0.8676 



The unit transport costs have been modeled here by means of a quadratic function 
shaping the dependency of the costs on loss factor for different route lengths 

c T = a- f 2 +b ■ f + c [ €/t ] . (3.10) 



The set of parameters ( a,b,c ) adopted for each distance value (/_) are in Table 3.3. 
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Table 3.3 Parameters values of Eq. (3.10) 



L 


a 


b 


c 


[km] 


re/ti 


re/ti 


re/ti 


200 


391.37 


-402.35 


108.69 


300 


370.86 


-354.62 


108.34 


400 


350.35 


-306.9 


107.99 


500 


329.84 


-259.17 


107.63 


1000 


227.28 


-20.54 


105.87 


2000 


22,168 


456.72 


102.35 


3000 


-182.94 


933.98 


98.82 


4000 


-388.06 


1411.2 


95.3 


5000 


-593.17 


1888.5 


91.77 


10,000 


-1618.7 


4274.8 


74.15 



By adopting Eqs. (3.8), (3.9) and (3.10), the logistics cost function can be 
reformulated as: 



— [a-f 2 +b-f + c\+ 



F,.=- 



Cu 



2 ■ m ■ L 



p H \t v f 2 + k 2 -f + k,\ 



f ■'Le r y i +f-'Le i -a i +'Le i -P i 



(3.11) 



As an example, the values of F L and of its three terms (transport, inventory, and 
environmental costs, respectively) against the loss factor value, for given values of 
p (0.5), L (200 [km]) and c/, (5000 [€/t year]) are plotted in Fig. 3.2. As it can be 
observed, it is possible to identify a proper / value allowing minimizing both 
logistics and environmental cost. 



L = 200 [km]; p = 0.5; * = 0.5; c h = 5000 [€/t] 




Fig. 3.2 Transport, environmental, holding costs and logistic cost factor Fi 
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In case of p = 0.5 and C/, = 5000 [€/t year], problem (3.6) has been solved for 
different L values (200^100-1000-2000 [km]). Results are in Fig. 3.3. As it can be 
observed, the optimal loss factor value (the loss factor which minimizes the global 
logistic function) decreases from about 0.6 to less than 0.1 with the increase of the 
route lengths from 200 [km] to over 1000 [km]. It can be concluded that in case of 
long route lengths, slow transport means, characterized by an average speed less 
than 50 [km/h], reveal a sustainable choice. In case of L > 500 [km], different 
transport means (rail: /= 0.035, or ships: /= 0.010) are characterized by 
approximately the same global logistic costs. On the contrary, in case of short 
distances ( L < 500 [km]), a proper means of transport (tir < 7.5 [t]:/= 0.475, 
small tir < 3.5 [t]: /= 0.718) allows minimizing problem (3.11). Finally, it can be 
observed how transport by airplane (f= 0.989) does not represent a sustainable 
choice for any distance considered, since F L assumes the highest value for any 
route length. 



p = 0.5; k = 0.5; c h = 5000 [€/t] 



0.350 • 
0.300 • 
0.250 - 
0.200 - 
0.150 - 
0.100 ■ 
0.050 ■ 
0.000 - 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

| 200 [Km] 400 [Km] 1000 [Km] 2000 [Km]| 

Fig. 3.3 F l values for different route lengths and / values 



As discussed above, the solution of problem (3.6) allows identifying the 
optimal loss factor value {[opt) and hence the optimal transport means for a given 
transport distance L. Each transport means is characterized by a proper average 
transport speed v*(f OP1 ). 

By introducing v*(f 0 pT) in Eq. (3.3), the Sustainable Order Quantity, SOQ, 
minimizing the overall logistic cost function (Eq. 3.1), is obtained as: 



p ■ H -V* (f OPT ) 

Results of a sensitive analysis carried on SOQ/G ratio and /opt values varying 
the transport distance L, for different p values (H = 3520 [h/year]; k = 0.5; 
Ch = 5000 [€/t]) are in Fig. 3.4. (SOQ/G) 1 ratio represents the number of annual 
orders. 
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(a) SOQ/G 





ap = 0.25 up = 0.5 ap = 0.75 sp = 1.0 



L [km] 



Fig. 3.4 SOQ/G (a) and j, opt (b) versus transport distance L for different p values 

It can be observed how, according to a sustainable logistics policy, in case of 
transport distance of more than 500 [km], slow transport means (low /opt values) 
and high order quantity should be adopted; in case of route lengths shortest than 
500 [km], the solution move towards small SOQ/G values (see Fig. 3.4a) with 
medium average speed transport means (e.g. light trucks) (see Fig. 3.4b). 

Parameter p (see Fig. 3.) represents the ratio between the supply lead time (LT) 
and the consumption time ( CT ). It proves to slightly affect f OPT values (see Fig. 
3.4a): even in the case of small supply lead time (p = 0.25), the proposed model 
does not identify high speed transport means as a sustainable choice because of 
the high environmental costs. On the contrary, SOQ/G value increases with the 
decreasing of the allowed (maximum) supply lead time/consuming time ratio. 
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Model Application 

For given values of the route length L and of the holding cost c h , problem (3.6) 
can be solved for different p values. The model provides the optimal solution of 
the logistics problem that is the couple ( SOQ ; transport means) minimizing the total 
logistics costs. As an example, in Fig. 3.5a solutions obtained for different p 
values in case of a transport distance L = 400 [km] and of a unit holding cost 
value c/, = 5000 [€/t-year] are showed. For each p value considered, the model 
identifies in trucks with different capacity (C) the optimal solution. In Fig. 3.5a, 
the corresponding SOQ values and the capacity of the truck (C) identified as 
optimal transport means are summarized. In case of L = 440 [km] the optimal 
solution consists of a truck with C < 7.5 [t], in the optimal order quantity 
SOQ = 0.0041 -G, and in a p value of 0.5, that is in a reorder quantity r = 0.5- SOQ. 

On the contrary, in case of L = 1000 [km], the optimal means of transport 
identified for any p value is a ship (f= 0.010); however, different SOQ values are 
obtained varying parameter p (see Fig. 3.5b). 

(a) F l L = 400 [km]; c„ = 5000 [€/t year] 




( b ) f l L = 1000 [km]; c„ = 5000 [€/t year] SOQ/G 




Fig. 3.5 Fl, SOQ/G, and/bir values for c/, = 5000 [€/t-year] in case of (a) short distances 
(L = 400 [km]) and (b) long distances ( L = 1000 [km]) 
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In order to assess the feasibility of the theoretical solution identified by the 
model, for each p value a final check is required: the allowed supply lead time 
[LT all = L/v(f OPT )], has to be greater than the average supply lead time (LT), being 
LT evaluated on the basis of the average speed of the transport means selected 
( LT = L / v). 

As an example, in case of p = 0.5, the limit (allowed) values of the supply lead 
time are listed in Table 3.4 for different route lengths (k = 0.5; C/, = 5000 
[€/t-year]). 

Table 3.4 Allowable lead time ( LT all ) values obtained by the model in case ofp = 0.5 ( k = 0.5; 

Ch = 5000 [€/t-year]) 



L [km] 


200 


300 


400 


500 


1000 


2000 


LT all [hi 2.5 


4.4 


6.6 


8.1 


16.2 


32.4 



In case of no transport means available satisfies the above constraint, three 
strategies could be adopted. 

(a) Adoption of a different reorder level, that is of a different sub-optimal 
solution; generally, this entails an increase in the logistic costs as per Fig. 3.5a; in 
case of long transport distance (L > 500 [km], see Fig. 3.5b), however, negligible 
differences are observed; as an example, in case of L = 400 [km], the sub-optimal 
solution 4 satisfies the lead time constraint ( LT all < LT) with an increase of less 
than 1 % in the logistic costs, as per data in Table 3.4. 

(b) In case no sub-optimal solutions are available, a faster available means of 
transport has to be selected with the result of an increase in the logistic costs; as an 
example, in case of Z, = 300 [km], it is not possible to find a feasible solution (see 
Table 3.5). For p = 0.5, by adopting a means of transport that is faster than the one 
identified by the model (solution 3mod. in Table 3.6) it is possible to satisfy the 
time constraint ( LT all < LT), but this leads to an increase in the total logistic costs 
of about 30 %. 

(c) In the case that there is no faster means of transport available, in order to 
satisfy the time constraint a lower k value is required, that is material handling, 
order management, and the quality control time have to be reduced. Generally, 
low decreases in k are possible. 
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Table 3.5 Solutions of the logistics problem (36) in case of L = 400 [km] (A: = 0.5; c* = 5000) 



Solution 


P 


SOO/G 


LT all 

Thl 


/(DPT 


Transport 

means 


LT 


~ f l 


1 


0.10 


0.0100 


3.54 


0.613 


Truck; 

C < 3 . 5 [t] 


6.00 


0.0438 


2 


0.25 


0.0061 


5.38 


0.530 


Truck; 

C < 7.5[t] 


7.50 


0.0415 


3 


0.50 


0.0041 


7.29 


0.476 


Truck; 

C < 7.5[t] 


7.50 


0.0411 


4 


0.75 


0.0033 


8.62 


0.449 


Truck; 

C < 7.5[t] 


7.50 


0.0412 


5 


1.00 


0.0027 


9.66 


0.430 


Truck; 

C < 7.5 m 


7.50 


0.0415 



Table 3.6 Solutions of the logistics problem (36) in case ofZ = 300 [km] (k = 0.5; c/, = 5000 

[€/t-year]) 



solution 


P 


SOQ/G 


LTall 

[h] 


foPT 


transport 

means 


LT 

[h] 


~ f l 


1 


0.10 


0.0078 


2.74 


0.606 


truck; 

C < 3.5[t] 


4.50 


0.0293 


2 


0.25 


0.0043 


3.81 


0.540 


truck; 

C < 7 . 5 [t] 


5.63 


0.0280 


3 


0.50 


0.0025 


4.48 


0.502 


truck; 

C < 7.5[t] 


4.48 


0.0277 


3mod. 


0.50 


0.0025 


4.48 


0. 7176 * 


truck ; 

C < 3.5[t] 


4.48 


0.0357 


4 


0.75 


0.0020 


5.21 


0.485 


truck; 

C < 7.5[t] 


5.63 


0.0278 


5 


1.00 


0.0016 


5.54 


0.474 


truck; 

C < 7.5[t] 


5.63 


0.0279 
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3.1.2 Environmental and Social Costs 

The main goal of the model application proposed in this Section is to evaluate how 
different internalization strategies can promote the shift from the adoption of 
transport means generating high external costs toward transport means 
characterized by better environmental and social performances. 

Starting from Eq. (3.1) and Eq. (3.3) and explicitly considering the order costs, 
the annual logistics cost function can be evaluated as: 






/r T \ 1 C/i ’l* 

m-c T (fL )-\ — 



2 p-H-v 



+ C r 



pH -v 
L 



(p 

^ r 



(3.13) 



Table 3.7 summarizes the unit external costs adopted in this Section: cost data 
for water, rail and road freight transport come from [3], cost data for air freight 
transport comes from [4]. Unit external costs of road transport are provided for 
both LDV (Light Duty Vehicles) category, with a GWR (Gross Weight Rate) less 
than 3.5 [t], and HDV (Heavy Duty Vehicles) category, with a GWR greater than 
3.5 [t]. The LCA category considers the up- and down processes external costs. 
Cost data have been updated to [€ 2 oi 3 /fkm] values by means of discounting 
indexes available in [5]. 

Table 3.7 Unit external costs [€ 2 oi 3 /fkm] of different transport means 





Water- 

borne 


Rail 

electric 


Rail 

diesel 


LDV 


HDV 


Airplane 


Cost category 






[€ 20 i3/fkm] 






Accidents 


0.0 


0.2 


0.2 


62.3 


11.3 


0.0 


Air Pollution 


6.0 


1.0 


1.9 


19.8 


7.4 


20.7 


Noise 


0.0 


1.1 


1.1 


7.0 


2.0 


11.8 


Congestion 


0.0 


0.0 


0.0 


107.0 


28.5 


0.0 


GW 


4.0 


0.0 


4.3 


49.3 


10.9 


325.9 


LCA 

Other external 


1.4 


4.4 


5.7 


15.8 


3.3 


10.2 


cost 


1.0 


0.6 


0.6 


7.1 


2.8 


5.0 


Overall 


12.4 


7.3 


13.7 


268.3 


66.2 


373.7 



By adopting average values of the loss factor for the different transport 
modalities in Table 3.7, it is possible to shape the dependency of the overall unit 
external costs on the loss factor by means of a quadratic function: 



£ ex =8-f 2 +p-f. 



( 3 . 14 ) 
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Being: 



^EX ~ £ EX (■/") ' m 'G, 



(3.15) 



by adopting quadratic functions to shape the dependency of the speed (see Eq. 
3.8) and of the unit transport costs (see Eq. 3.10) on the loss factor as in the 
previous Section, the logistics cost factor F L is: 



• m 



■\_a-f 2 + b-f + c~\ + — 



f l =- 



P-H [k r f 2 +k 2 -f + k 3 ] 



+ 



2 'r° P f {^/ 2+ ^/ + ^] + — {^/ 2 + A 7 ] 

\j * Ctj * -L/ C, 



2 ■ m ■ L 



>.( 3 . 16 ) 



A Numerical Example 

In this case study, updated data on loss factor (/) values as well as on transport 
speed (v) have been adopted. Data are in Table 3.8. 



Table 3.8 Loss factor (/) and average speed of transport (v) values for different means of 
transport; [DB1] = [1], 2012; [DB2] = [6], 2012; [DB3] = [7], 2007; [DB4] = [8], 2012 



Transport means 


Database 


/ 


Database 


V 

[km/h] 


Ship — oversee 


[DB1] 


0.010 


[DB3] 


24.7 


Ship — tank 


[DB1] 


0.021 


[DB3] 


18.5 


Ship — coast 


[DB1] 


0.033 


[DB3] 


17.3 


Rail — diesel 


[DB1] 


0.033 


Assumed 


46.7 


Rail — electric 


[DB1] 


0.037 


Assumed 


46.7 


Ship — inland 


[DB1] 


0.048 


[DB3] 


24.7 


Truck 34 - 40 [t] 


[DB1], [DB2] 


0.063 


[DB4] 


56.5 


Truck 28 - 34 [t] 


[DB1], [DB2] 


0.069 


[DB4] 


56.5 


Truck 26 - 28 [t] 


[DB1], [DB2] 


0.119 


[DB4] 


56.5 


Truck 20 - 26 [t] 


[DB1], [DB2] 


0.136 


[DB4] 


56.5 


Truck 14-20 [t] 


[DB11, [DB2) 


0.127 


[DB4] 


64.0 


Truck 12 — 14 [t] 


[DB1], [DB2] 


0.139 


[DB4] 


64.0 


Truck 7.5 - 12 [t] 


[DB1], [DB2] 


0.149 


[DB4] 


64.0 


Truck 3.5 - 7.5 [t] 


[DB1], [DB2] 


0.165 


[DB4] 


64.0 


Truck 3.5 - 7.5 [t] a 


[DB1], [DB2] 


0.156 


[DB4] 


64.0 


Truck <3.5 [t] b 


[DB1], [DB2] 


0.335 


[DB4] 


72.7 


Truck <3.5 [t] c 


[DB1] 


0.554 


[DB4] 


72.7 


Airplane 


|DBI| 


0.984 


Assumed 


466.7 



“Truck operated without append 

b Light good vehicles 

c Vans 
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Data in [8] on free flow vehicle speeds observed in UK in 2011 have been 
elaborated in order to obtain an average speed value for each truck category (see 
Table 3.9). The actual speed value (v ACT ) for each truck category has been 
obtained from observed data by weighing the values observed on the different 
road types (motorways, dual carriageways, single carriageways) with the 
corresponding observations number. The related average value of the transport 
speed has been obtained using a k factor of 0.5. 

Data on ships cruise speeds as in [7] have been adopted. For railroads and 
aircraft transport means, a cruise speed of 70 [km/h] and 700 [km/h] have been 
assumed, respectively. 
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By relating the loss factor values of different transport means with their 
characteristic average speed in transport (obtained assuming k = 0.5), the data set 
in Table 3.8 has been obtained. Results of a regression analysis carried out on the 
data set showed that the relationship between the speed of transport and the loss 
factor value of the transport means adopted is best shaped by a quadratic function: 

v = k i -f 2 +k 2 -f + k 3 , (3.17) 

with k, = 557.6; k 2 = -150.4; k 3 = 52.7 (R 2 = 0.935). 

Starting from the data set in Table 3.8, average values of the loss factor for the 
different transport modalities have been evaluated. Values adopted are in Table 
3.10. 



Table 3.10 Average values of the loss factor for different transport modalities 



Transport means 


f 


Sea Ship 


0.021 


Diesel Train 


0.033 


Electric Train 


0.037 


Inland Ship 


0.048 


HDV 


0.125 


LDV 


0.444 


Airplane 


0.984 



A regression analysis carried out on data in Table 3.7 led to evaluate 
parameters values of Eq. (3.14) as: 8= -607.17, // = 979.22 (R 2 = 0.9718). 

Table 3.11 summarizes the data set adopted for the numerical experiment. 

In order to evaluate the optimal transport means and the optimal lot-size value 
minimizing the logistics function F L , problem (3.6), considering also Eq. (3.3), has 
been solved. For each transport distance considered both the solution obtained in 
case of internalization of the whole external costs (F LSUST in Table 3.12, SOQ SU st 
in Table 3.13 and r SUST in Table 3.14) and the solution obtained in case of no 
internalization of external costs (F IEC0N in Table 3.12, SOQecon in the Table 
3. 13 and r EC0N in Table 3.14) are showed. In brackets, the optimal transport means 
is specified (IV = Inland Vessel, ET = Electric Train, SS = Sea Ship). 

The optimal lot size and the reorder level, according to the model defined, 
proved to be independent on the external costs. For each transport means, SOQ 
and r values depend only on the transport distance considered (see Table 3.15 and 
Table 3.16). 
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Table 3.11: Data set adopted for the numerical experiment 



Parameters 


Values 


Measurement units 


G 


9224 


[unit/year] 


H 


3520 


[h/year] 


m 


0.5 

200; 300; 400; 500; 1000; 


[kg] 


L 


2000; 3000; 4000; 5000; 
10,000 


[km] 


c H 


4.65 

See Eq. (3.10) and Table 


[€/unit-year] 


c t 


3.3 


[€/t] 


Co 


1 


[€/order] 



Sensitivity Analysis 

In order to evaluate the impact of different internalization strategies on specific 
logistics costs [€/kg-km], the problem (3.6) has been solved for different extra 
costs charged in the logistics costs function. Ten different internalization strategies 
have been considered, with a percentage of charged extra costs ranging between 
0 % and 200 % of the whole external costs (100 %-case values are in Table 3.7). 
The specific logistics cost is defined here as: 



<t>L 



G-m-L 



(3.18) 



Starting from solutions obtained by the model, the increase of the specific logistics 
cost (compared with the base case of zero external costs charged) have been 
computed for each case. Values are in Table 3.17, where in brackets there is the 
optimal transport means identified by the model in each case. 



Table 3.12 Fl.econ and F^susr for different transport distances 
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Table 3.15 Sustainable Order Quantity (SOQ) for different transport means and distances 
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For the shortest transport distance considered (200 [km]) road transport is the 
optimal solution for each of the internalization strategies considered. An external 
costs charging, in the logistics function, of more than 160 % of the whole external 
cost generated from the transport causes a shift from LDV to HDV transport 
modality. 

In case of short distances ( L = 300, 400, and 500 [km]), road transport means 
are the optimal solution of the problem not for all the internalization strategies 
considered. Depending on the transport distance, the increase of the percentage of 
the whole external costs charged in the logistics cost function promotes the shift 
toward more sustainable transport modes (inland vessel, electric train). 

Road transport is not solution of the problem (3.6) in case of transport distances 
longer than 1000 [km] for all the internalization strategies considered. 

Fig. 3.6 shows the trend of the specific logistics costs for all the transport 
means considered in case of four different distances ( L =200, 500, 1000, 2000 
[km]) and three different internalization strategies (0 %, 100 %, and 200 %). In 
Fig. 3.7 the specific logistics costs for different transport distance (L) and different 
internalization policies adopted are depicted. 




Sea Ship Diesel Electric Inland HDV LDV Airplane 

Train Train Vessel 



Vi [€/t km] (l=soo 

1,6 
1,4 
1,2 
1,0 
0,8 
0,6 
0,4 
0,2 
0,0 

Sea Ship Diesel Electric Inland HDV LDV Airplane 
Train Train Vessel 
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<Pi [€/t km] 




■ 100 % 
-200% 



i Diesel Electric Inland HDV 
Train Train Vessel 




— 0% 

• 100 % 
—200% 



Fig. 3.6 Specific logistics cost for different transport means and different internalization 

strategies. 
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Fig. 3.7 Specific logistics cost percentage increase compared to the economic case for 
different transport distance ( L ) and two different internalization strategies including all the 

external costs categories 



As shown in Fig. 3.7, for a given distance, the adoption of a more severe 
internalization strategy results in an increase of the specific logistics cost which 
significantly differs depending on the transport means considered. As an example, 
for a transport distance of 1000 [km] and with an increase from 100 % to 200 % of 
the percentage of external costs charged, the increase in the specific logistics cost 
is below of the 8 % for the sea ship and electric train transport mode. On the 
contrary, substantial differences have been observed in case of road transport 
modalities: 28 % in case of HDV and 57 % in case of LDV. 

GW and LCA costs assume the larger values in the economic evaluation of the 
external costs of transport: the sum of these two categories weights up to 90 % (air 
transport) of the total external costs. Moreover, in the evaluation of sustainability 
in the transport sector, often only the carbon emissions are considered. A 
sensitivity analysis has been carried out in order to evaluate the effects of different 
internalization strategies if only GW and LCA external costs are charged in the 
logistics cost function. Results obtained are in Table 3.18 and in Fig. 3.8. 



3.1 Deterministic Demand and Lead Time 65 



< 

U 

X 

-a 

u 



o 

13 

O /■ 



C3 , — i 

<L> 

O > 

«j 2, 
QJ 

■a 

0> --H 

•a OX) 
& Vh 
J3 ctf 
c« *d 
C ° 
O i2 



t3 x 

• - pq 



o o 

(D g 
c/a O 
c3 O 



‘-d o3 

’W) £ 

O o 



<L> o 
Oh t>p 

03 1) 



& 8 

cS 

a 2 
a £ 

• S <L> 



Oh 

00 



-D 

C3 

H 





































C/T 


CO 


CO 

GO 


S? 

GO 


co' 


co' 


go' 


go' 


go' 


go' 


go' 






o 


GO 


GO 


GO 


GO 


GO 


GO 


GO 






o 

o 


w 

N® 


w 

X® 


X® 


X® 


^P 

o x 


^P 

o x 


© x 


^P 

© x 


^P 

o x 


^P 

o x 


^P 

© x 






o' 


cN 


o'- 


o'- 


0 s - 


O 


CO 


in 


00 


O 


CO 


in 






’ 1 


o 


CO 


in 


00 


* 1 


* 1 


’ 1 


’ 1 


<N 


<N 


CN 




























/ s 








s 


s 


,. — s 


,. — s 


^ N 


GO 


GO 


GO 


GO 


GO 


GO 








GO 


GO 


GO 


GO 


GO 




GO 






GO 


GO 














GO 




w 




w 


w 




s — f 






o 


w 


w 


w 




w 


















o 












o s 


© s 


© s 


© s 


© s 








o 












O 


(N 


^t 


\D 


OO 


O 






IT) 


O 


<N 


T ^t 


NO 


00 


’ 1 


’ 1 


* 1 


' 1 


’ 1 


(N 




















^ ^ 


^ ^ 





^ ^ 


/ s 








s 


,. — s 


,. — s 


,, — s 


^ N 


^ N 


GO 


GO 


GO 


GO 


GO 








GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 








GO 


GO 


GO 




GO 


GO 










s — ** 






o 








w 




















o 

o 


N® 

© x 


N® 

o x 


N® 

o x 


p" 

o x 


o x 


o x 


© s 


© s 

CO 


© s 


© s 

\D 


00 








O 


<N 


T ^t 


•n 


r- 


Os 


’ 1 


’ 1 




' 1 


r— 1 




























/ s 








— s 


^ s 


,. — s 


^ N 


^ N 


^ N 


^ N 


GO 


GO 


GO 


GO 








GO 


GO 


GO 


GO 


GO 


GO 


GO 






GO 


GO 








GO 






GO 








w 


w 




s — f 






o 

o 




w 


w 




w 


w 


w 


N® 

o x 


vP 

o x 


vP 

o x 


P 

o x 






o 


















(N 




«n 






CO 


O 


(N 


CO 


«n 


NO 


oo 


On 


’ 1 


* 1 


’ 1 


’ — 1 


























^ ^ 




























GO 


GO 








00 


GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 


GO 








GO 


GO 






GO 


GO 


GO 


GO 


GO 




s — f 






o 

o 






w 


w 












-P 

o x 


P 

o x 






o 






















<N 






<N 


o 


,_H 


<N 




«n 


NO 


r- 


00 


On 


1—1 


' 


1 — 1 




























£ 










/"x 


f — ' 














/ — N 


-x 








H 


H 


H 


H 


H 


H 


H 


H 


H 


H 










PP 


PP 


PP 


W 


W 


PP 


PP 


W 


PP 


W 






o 




























o 

o 


p 

o x 


? 


x© 

o x 


nP 

o x 


NP 

o x 


o x 


N® 

© X 


N® 

© X 


N® 

o x 


N® 

o x 


P 

© x 








O 


1—1 


1 


<N 


<N 


CO 


CO 






«n 


«n 





































> 

Q 

hP 


> 

Q 

_p 


> 


> 


> 


> 


> 


> 


> 


> 


> 








Q 


Q 


Q 


Q 


Q 


Q 


Q 


Q 


Q 












hP 


X 


X 


X 


X 


X 


X 






o 








sP 

o x 


^P 

o x 


? 


s? 


s? 


£ 


5? 


£ 






o 








NO 


<N 


in 


VO 


t~ 


00 


ON 








IT) 


o 


in 


* 1 


’ 1 


<N 


<N 


<N 


<N 






m 












^ ^ 


^ 


^ 





^ ^ 


^ ^ 


^ ^ 


, , 










> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 








Q 


Q 


Q 


Q 


Q 


Q 


Q 


Q 


Q 








Q 


Q 








X 


X 


X 


X 


X 


X 








|_P 


|_P 
























o 






X® 

o x 


o x 


sP 

o x 


s= 

o x 


v? 

o x 


s= 

o x 


o x 


© x 


o x 






o 






o 


in 


O 


«n 


O 


«n 


O 


CO 










o 


<n 


* 1 


’ 1 


fN 


(N 


CO 


CO 






-St - 





































> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 








Q 


Q 


a 


Q 


Q 


Q 


Q 


Q 








Q 


Q 


Q 


►J 


hP 


X 


X 


X 


X 


X 


X 








hP 


hP 


hP 




































'P 






s= 


s= 






o 








© x 


o x 


© x 


© x 


© x 


© x 


© x 


© x 






o 








CO 


00 


(N 


r- 


i— i 


N£3 


O 


in 






CO 


o 




On 


’ 1 


’ 1 


<N 


(N 


CO 


CO 


^J- 


^t 





































> 


> 


> 


> 


> 


> 


> 


> 


> 


> 


> 








Q 


Q 


Q 


Q 


Q 


Q 


Q 


Q 








Q 


Q 


Q 






X 


X 


X 


X 


X 


X 








hP 


hP 


hP 










































n? 


N? 






o 








o x 


o x 


© x 


© x 


© x 


© x 


© x 


© x 






o 








i— i 


«n 


00 


(N 


N£3 


On 


CO 


r- 






CN 


o 






,_H 




1—1 


<N 


(N 


(N 


CO 


co 


OX) 


























a 














N® 


nP 


N® 


N® 


v® 


P 








p 

o x 


X© 

o x 


nP 

o x 


nP 

o x 


0 s 

O 


0 s 

O 


0 s 

O 


0 s 

O 


0 s 

O 


0 s 

O 


CO 


►> 

JD 




O 


O 


O 


O 


O 


<N 




NO) 


00 


O 


U 


o 


<N 




NO 


00 


* 1 


’ 1 


’ 1 


* 1 


* 1 


<N 



66 3 SOQ Model Formulations 



<p L (Sea Ship) 



> t (Electric Train) 



50 

' 200% 40 



’--200% 
► 100 % 



o ^ ' 

^ ^ ^ 

L [km] 

% increase <p L (HDV) 



1 

■P 



% increase 
80 



v v ^ <r 
Z. [km] 

<Pl (LDV) 



- 200 % 
► 100 % 



- 200 % 
► 100 % 



■p # # •? ^ 



L [km] 






4 # & <? / # # # ^ 

L [km] 



Fig. 3.8 Specific logistics cost percentage increase compared to the economic case for 
different transport distance ( L ) and two different internalization startegies charging only 
GW and LCA external costs categories. 



The trends in Fig. 3.8 are similar to those of Fig. 3.7 but with a smaller specific 
logistics costs increase (8 % in case of FID V, 24 % in case of LDV and below the 
4 % in case of electric train and the sea ship). 

The internalization strategies charging only the costs of GW and LCA of 
transport reveal to be ineffective in terms of more environmental and social 
friendly transport means selection in case of short transport distances (L < 500 
[km], see Table 3.18). 



3.2 Stochastic SOQ Model 

The formulation of the SOQ model proposed in Sect. 3.1 allows identifying the 
optimal lot-size and the transport means minimizing logistic and environmental 
costs of transport in case of a deterministic product demand. In case of stochastic 
variability of product demand and/or supply lead time, shortage costs, as well as 
safety stock costs, have to be considered in the logistics cost function (see Eq. 
3.19). 



+ [€/year] 



(3.19) 
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Usually, safety stock is sized searching for the trade-off between holding and 
shortage costs due to stock out events. In this Section, two extension of the SOQ 
model in case of stochastic variability of product demand and of supply lead time 
are discussed. Safety stock and shortage costs are considered in the logistics cost 
functions proposed. Analytical models provide order quantity, transport means 
and safety stock size allowing minimizing logistics and environmental costs in 
case of stochastic variability of product demand or of supply lead time. Results of 
the application of the models to a full scale case study from automotive industry 
are presented and discussed. 



3.2.1 Product Demand Uncertainty 

The following assumptions are assumed in defining the SOQ model in case of 
stochastic variability of the product demand: 

1. the expected annual demand ( G ) is known; 

2. the product demands in each time period, D„ are independent stochastic 
variables; 

3. in each z'-th period of the lead time ( LT ), product demand is characterized by 
the same expected value and by the same standard deviation: 



E(Dj) = E(D, ) = ... = E(D lt ) = E(D) = G/H, 



_2 _ _2 _ 



_ 2 _ 2 

v DLT °D ' 



(3.20) 

(3.21) 



As far as concern the transport speed (v) and the evaluation of the lead time 
(LT), the same assumptions as in Sect. 3.1.1 are adopted. 

In case of a product demand varying stochastically, and environmental costs of 
transport considered, the logistic cost function (3.19) can be written as: 



= c h 



2 Q + SS 



+ G ■ m ■ c T H 

Q 



■N s + G-m-L-Y j£r e(f). (3.22) 



In Eq. (3.22) the safety stock level, SS, consistent with an assigned service level, 
SL, is evaluated as: 

SS = LT-[D* (3.23) 

where: 

LTD * 

SL = prob | D tot <LT -D* j = J pdf(D TOT )dD TOT . 



-oo 



(3.24) 
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The number of stock out events is computed (see Sect. 2.4.4) by means of 

+00 

Ng — J pdf(D TOT ) • (D tot — Dj,q T ) • dD T Q T , (3.25) 

D T qt 



where 



pdf (D T ot ) 





(P TOT -LT ■ E (D)f 
2-LT-ol \ 



and 



D^ot — D\ T 4- .. . -F D • (3 .27) 

In Eqs. (3.23)-(3.25) D* and D TOT are the maximum demand value of the z'-th 
period and the maximum demand value on LT periods not causing a stock-out in 
that period and in LT periods, respectively. From Eq. (3.27) it is obtained: 

LTD* = Dj'qy . 

The optimal values of loss factor (/opt ), order quantity (SOQ), and the 
corresponding values of the reorder level (rtfopr)), and of the safety stock 
(SSpopp)) jointly minimizing the logistics and external costs can be obtained by 
solving problem (3.28) or the equivalent problem (3.29): 



min <P, 

SOOf 


(3.28) 


min F, 

P.f 


(3.29) 



By solving problem (3.27) or (3.28) the evaluation of the optimal loss factor 
value led to jointly identify the optimal order quantity (SOQ), and the optimal 
safety stock level (.S'.Sj as: 



SOQ = °' L (3.30) 

Popt ' n ' v \Jopt) 

SS = -^-\D\f OPT )-E(D)] (3.31) 

V \J OPT) 



Case Study 

The model has been applied to a case study from automotive industry [9, 10]. 
The case refers to a multi-site manufacturing system producing breaking 
equipment. The supply chain consists of three production sites (see Fig. 3.9): sites 
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1 and 2 are responsible of producing semi-finished products for sites 2 (PIC, P2C) 
and 3 (P1B, P2B); in site 3, three finite products (PI, P2, P3) are assembled for a 
single customer. Sites 2 and 3 produce additional products (PX, PY) required by 
external customers of the aftermarket starting from semi-finished product (PXf, 
PY f) externally supplied. The products demand of the main customer, as well as 
of the external customers, are uncertain and vary stochastically. 

The logistics problem (3.28) has been solved for PXf product in site 1. The PXf 
product demand is characterized by an expected value (E(D)) of 2.62 [unit/h] and 
a coefficient of variation cv of 0.26. The number of working hours ( H) in site 1 is 
3520 [h/year]. The maximum production capacity of product PX in site 1 is 18.40 
[unit/h]. This value is a theoretical limit for the production capacity, since it is 
obtained without considering the set-up time. 

The unit holding cost of product PXf in site 1 (c h ) is 15.49 [6/unit], and the 
extra cost generated in case of stock out (c s ) is 10.06 [€/unit] (c/c/, = 0.65). The 
mass (m) of the product PXf is 0.5 [kg] (Fig. 3.9). In evaluating (3.23) and (3.24), 
a lower bound value of 0 [unit/h] and an upper bound value of 18.40 [unit/h] have 
been considered, respectively. The analysis has been limited to road transport 
because of geographical constraints. Due to unavailability of data, order costs 
have not been considered. 




Fig. 3.9 The supply chain of a multi-site manufacturing system 

In this case study, updated data on emission factors (e,) have been adopted. 
Data are in Table 3.19. 

Loss factor (/) and average transport speed (v) values adopted are those of Sect. 
3.1.2. As in Sect. 3.1.2, average speed of transport has been shaped as a quadratic 
function of the loss factor: 



v = k 1 -f 2 +k 2 -f + k 3 



(3.32) 



with k i = 557.6; k 2 = -150.4; k 3 = 52.1. 

Equivalent emission values in Table 3.19 have been obtained starting from data 
on emission factors of different transport means in [1, 6]. Data allowed evaluating 
the impacts on the environment in terms of Global Warming (GW), Acidification 
(AP) and Tropospheric Ozone Precursor Potential (TOPP) of different means of 
transport. Starting from the unit emission [kg/t-km] of different pollutants, the 
impacts have been measured by evaluating the equivalent emissions of CO 2 (GW), 
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S0 2 (AP), and NMVOC (TOPP). In case of GW and AP, classification factors as 
in [11] have been adopted; in case of TOPP classification factors as in [12] have 
been adopted (see Table 3.20). 



Table 3.19 e ; values for different means of transport; [DB1] = [1], 2012; [DB2] = [6], 2012 



Transport means 


f 


Database 


[kg/t-km] 


(SCMen 

[kg/t-km] 


NMVOCen 

[kg/t-km] 


Ship — oversee 


0.010 


[DB1] 


0.027 


0.000 


0.000 


Ship — tank 


0.021 


[DB1] 


0.018 


0.001 


0.000 


Ship — coast 


0.033 


[DB1] 


0.044 


0.000 


0.000 


Rail — diesel 


0.033 


[DB1] 


0.032 


0.000 


0.000 


Rail — electric 


0.037 


[DB1] 


0.083 


0.001 


0.000 


Ship — inland 


0.048 


[DB1] 


0.044 


0.000 


0.001 


Truck 34-40 |t| 


0.063 


[DB1], [DB2] 


0.054 


0.000 


0.000 


Truck 28 - 34 [t] 


0.069 


[DB1], [DB2] 


0.059 


0.000 


0.001 


Truck 26 - 28 [t] 


0.119 


[DB1], [DB2] 


0.101 


0.000 


0.001 


Truck 20 - 26 [t] 


0.136 


[DB1], [DB2] 


0.116 


0.001 


0.000 


Truck 14-20 [t] 


0.127 


[DB1], [DB2] 


0.109 


0.001 


0.001 


Truck 12-14 [t] 


0.139 


[DB1], [DB2] 


0.119 


0.001 


0.001 


Truck 7.5 — 12 [t] 


0.149 


[DB1], [DB2] 


0.127 


0.001 


0.001 


Truck 3.5 - 7.5 [t] 


0.165 


[DB1], [DB2] 


0.142 


0.001 


0.002 


Truck 3.5 - 7.5 [t] a 


0.156 


[DB1], [DB2] 


0.136 


0.000 


0.001 


Truck <3.5 [t] b 


0.335 


[DB1], [DB2] 


0.291 


0.001 


0.001 


Truck <3.5 [t] c 


0.554 


[DB1] 


0.442 


0.003 


0.005 


Airplane 


0.984 


IDB11 


0.828 


0.005 


0.005 



“Truck operated without append 

b Light good vehicles 

“Vans 



Table 3.20 Classification factors adopted 



Burdens 




Impact categories 


GW 


AP 


TOPP 


ch 4 


21 




0.014 


CO 






0.11 


co 2 


1 






H 2 S 




1 




HC1 




0.88 




HF 




1.6 




n 2 o 


320 






NH, 




1.88 




NMVOC 






1 


NO x 




0.7 


1.22 


so 2 




1 
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The regression analysis on equivalent emissions data in Table 3.19 led to model 
the three unit mass of equivalent emissions e, [kg/t-km] as a linear function of the 
loss factor: 



e=a,-f + P i i = (C0 2 ) eq , (S0 2 ) eq , (NMVOC) eq . . (3.33) 



Regression parameters a, and ft, and corresponding R 2 values, as well as the 
monetary cost per unit mass emission of the z'-th pollutant c, as in [13] are in table 



3.21. 



Table 3.21 Regression parameters and unit monetary costs for the impact category considered 



Impact 


a, 


ft, 


R 2 


£i 


i 


[kg/t-km] 


[kg/t-km] 


H 


re/kgi 


GW 










(co 2 ) eq . 

AP 


0.8192 


0.0094 


0.9954 


0.05 


(S0 2 ) eq . 

TOPP 


0.0045 


0.0001 


0.9193 


4.00 


NMVOC ea . 


0.0051 


0.0002 


0.8443 


30.00 



Starting from data provided by an Italian logistic company (see Table 3.22) 
transport costs have been modeled here as: 

0 T = ^'[ z i ' e ^ Z2 m-Q + T ] -e^ Tl ■ L-m-o\ (3.34) 

In Eq. (3.33), the transport cost for a single order ( <p T , the terms in square 
bracket in Eq. 3.33) is obtained as the sum of the cost due to the number and the 
particular type of the Cargo Transport Unit (CTU), <p T ,cru > and the cost related to 
the covered distance and the mass of products transported, <p T y. 

<Pt = <Pt,ctu + Vtm (3.35) 

Both contributions have been shaped as a function of the loss factor value (/). 
In case of road transport, the CTU can be identified as the payload of the truck 
type considered. Results of a regression analysis carried out on the data provided 
by an Italian logistic company led to shape both contributions through exponential 
relationships (see Table 3.23): 



/' z 2 



Vt^tu ~ z i ' e 



■ m ■ Q 



(3.36) 
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( Pr,k= T \ ' ef Tl L m-Q 



(3.37) 



Table 3.22 Transport cost data adopted 



Transport 

means 


GVWR 

category 


/ 


C7T//Payl 

oad 


( Pt.ctu 


0T.CTU 


(Pr.k 


m 


re/CTU] 


re/ti 


[€/t-km] 


Ship 










107.8 


0.0032 


Rail 










88.0 


0.0114 


Truck 


34 - 40 [t] 


0.0633 


24.0 


975 


40.6 


0.0219 


Truck 


26 - 28 [t] 


0.1193 


12.0 


731 


60.9 


0.0328 


Truck 


7.5 - 12 [t] 


0.1488 


5.0 


488 


97.5 


0.0525 


Truck 


3.5 - 7.5 [t] 


0.1651 


3.0 


423 


140.8 


0.0758 


Airplane 










40.3 


0.2679 



Table 3.23 Regression parameters values of transport costs functions 



Transport means 


z; 




T, 




[€/t] 


Z2 


[€/t-km] 


t 2 


Ship 


107.800 


0.000 


0.032 


0.000 


Rail 


88.000 


0.000 


0.011 


0.000 


Truck 


17.780 


11.740 


0.010 


11.740 


Airplane 


40.300 


0.000 


0.268 


0.000 



In case of maritime, rail and airfreight transport, data from [14] have been 
adopted (see Table 3.22). Parameters values adopted in (3.36) and (3.37) are in 
Table 3.23. 

By adopting (3.25), (3.31), (3.32), and (3.34), the logistic cost function (3.19) 
can be rewritten as: 



0 , 



+LT ■ (D* - E(D) ) 



+ c. 



J_ G LT 

2 Hp 
pH 

'77^ Tr'J exp 

(LT) ' & D ‘ \ 27t d tot 



+ G 



■ m ■ [zj • i 



<f' z 2 



+ (T l -e fT *)-L\+ 



)-l\, 



(3.38) 



(P TOT - LT ■ E(D)) 2 
2-LT-al 



[D 



TOT D T ot \iDtoT + 



+ G-m-L- 



f ■'LE i -a i +Y.e i -p i 
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And the corresponding logistic cost factor as: 



F, = 



LT LT ■ (D* - E{D) ) 



Hp 



+ 2_.Cs 



2 G 
p-H 



+ 2 ■ c h ■ m ■ [zj ■ +(T x -e 



J-ti 



>4* 



(3.39) 



G c h (LT) V2 -o D -4l^ d‘„ 
2 



J exp 



(P TOT ~ LT ■ E(D)) 2 
2-LT-al 



1 Dmr A 



TOT ^ TOT r^TOT 



V 



■ 12-r/yj + 



+ m-L 

c h 



f ■ Z ■ O-i + Z 6/ ■ Pi 



The general logistic optimization problem (3.29) has been solved for different 
value of SL (90 %-95 %-99 %) and p (0.1 - 0.9) in order to evaluate the effect of 
the demand variability. Different values of the distance L (200-400-1000 [km]) 
have been considered. The upper limit value for the parameter p, due to 
organizational constraints, has been assumed equal to 0.9 instead of 1.0. Results 
obtained are in Table 3.24. 



Table 3.24 Results obtained in case of cs/c* — 0.65 







L = 


200 [km] 


sop 


SS 


SL 


foPT 


POPT 


F l 


[unit] 


[unit] 


0.90 


0.063 


0.4 


0.0161 


24.72 


1.67 


0.95 


0.063 


0.5 


0.0142 


19.77 


2.14 


0.99 


0.063 


0.9 


0.0122 


10.99 


3.04 






L = 


400 [km] 


SOQ 


SS 


SL 


foPT 


POPT 


F l 


[unit] 


[unit] 


0.90 


0.063 


0.6 


0.0191 


32.96 


2.35 


0.95 


0.063 


0.9 


0.0169 


21.97 


3.03 


0.99 


0.063 


0.9 


0.0152 


21.97 


4.29 






L = 


1000 [km] 


SOQ 


SS 


SL 


foPT 


POPT 


F l 


[unit] 


[unit] 


0.90 


0.063 


0.9 


0.0269 


54.93 


3.70 


0.95 


0.063 


0.9 


0.0252 


54.93 


4.77 


0.99 


0.063 


0.9 


0.0244 


54.93 


6.77 



The transport means with the highest capacity (the slowest ones) minimize 
logistical and environmental costs for each distance and service level considered 
even if the environmental costs are significantly lower than transport costs (10 %). 
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For each distance considered L, as expected, the safety stock size ( SS) increases 
with the increase of the service level (SL) fulfilled While for each SL value the SS 
size increases with L. 

For each distance considered, the service level (SL) differently influences the 
optimal order quantity. The optimal order quantity values (SOQ) decrease, with 
the increase of the service level, for low distances (L < 400 [km]); in fact, in case 
of high SL, smaller order quantities lead to an increase in shortage costs (see Eq. 
3.22) lower than the corresponding decrease of the holding costs. 

The increase of the service level achieved does not affect the solution of the 
problem in terms of order quantity for high distances (L = 1000 [km]). In this 
case, in fact, holding costs are significantly higher than shortage costs: the 
minimum of the logistic cost function is therefore obtained for the lowest values 
of the order quantity adoptable (when p = 0.9) compliant with lead time 
constraints. 

A sensitivity analysis has been carried out varying the unit shortage cost c s 
(c's = r • Ch, with v= 1, 2, 10) in order to evaluate the effect of the shortage costs 
on the solution of problem (3.28). 

The case of a unit shortage cost greater than the unit holding cost refers to low 
value-added products and high unit profit, as in case of mass production of high 
technology product. 

Results obtained in case of a distance equal to 200 [km] are in Table 3.25. 

In case of small transport distance, the solution of the problem is not affected 
by the unit shortage costs in terms of transport means selection since the optimal 
loss factor values does not change and is characteristic of the slowest and high- 
capacity means of transport. On the contrary, with the increase of the unit shortage 
costs, higher optimal order quantities are obtained as solution of the logistics 
problem. The increase of the transport distances (L = 400 [km], 1000 [km]) leads 
to the same remarks: the effect of increasing unit shortage cost on the solution of 
(3.29) is similar to the one observed in case of a short transport distances. In this 
case the increase in the optimal order quantity with the increase of unit shortage 
cost is less severe if compared with the short distance case due to the higher 
transport cost. 
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Table 3.25 Results of the c^/c* sensitivity analysis in case of L — 200 [km] 



SL = 90 % 











SOQ 


SS 


cJCh 


foPT 


POPT 


F l 


[unit] 


[unit] 


0.65 


0.063 


0.4 


0.0161 


24.72 


1.67 


1.00 


0.063 


0.3 


0.0175 


32.96 


1.67 


2.00 


0.063 


0.2 


0.0206 


49.43 


1.67 


10.00 


0.063 


0.1 


0.0336 


98.87 


1.67 






SL 


= 95 % 


SOQ 


SS 


cjc h 


foPT 


POPT 


F l 


[unit] 


[unit] 


0.65 


0.063 


0.5 


0.0142 


19.77 


2.14 


1.00 


0.063 


0.4 


0.0152 


24.72 


2.14 


2.00 


0.063 


0.3 


0.0172 


32.96 


2.14 


10.00 


0.063 


0.1 


0.0266 


98.87 


2.14 






SL 


= 99 % 


SOQ 


SS 


cjc h 


foPT 


POPT 


F l 


[unit] 


[unit] 


0.65 


0.063 


0.9 


0.0122 


10.99 


3.04 


1.00 


0.063 


0.9 


0.0125 


10.99 


3.04 


2.00 


0.063 


0.8 


0.0133 


12.36 


3.04 


10.00 


0.063 


0.3 


0.0168 


32.96 


3.04 



The effects of the demand stochastic variability on the solution of Eq. (3.29) 
have been evaluated through a sensitivity analysis. The logistics problem has been 
solved in case of a cv values of 0.5 and 0.8, cslch = 0.65 and L = 200 [km] (Fig. 
3.10 and Fig. 3.11). 

SOQ 




Fig. 3.10 SOQ vs. cv values in case of L = 200 [km] and cs/ci, = 0.65 
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ss 




Fig. 3.11 SS vs. cv values in case of L = 200 [km] and cs/ch = 0.65 



Again, for each service level ( SL ) considered and even in case high cv values, 
the slowest transport means minimize the logistic cost function. As expected, the 
increase in the demand variability leads to an increase of the SS level required. 
The demand variability also affects the optimal order quantity values. The increase 
in the optimal SOQ, however, reveals to be dependent to the service level 
fulfilled: in case of lower SL values, a greater increase of the SOQ values is 
observed. 



3.2.2 Lead Time Uncertainty 

The SOQ model is here defined in case of stochastic variability of supply lead 

time, under the following assumptions: 

a. the product demand ( D ) is deterministic and stationary, D = G/H; 

b. orders do not cross; 

c. backordering is not allowed; 

d. the expected value of the supply lead time ( E(LT )) is evaluated as the sum of 
the expected transport time ( E(T t )) and the expected value of the time 
required for the material handling, order management and quality control 
(E(T l ))\ E(T l ) is evaluated as a fraction ( k< 1) of E(T t ); consequently, speed 
of transport is obtained from the free flow or cruise speed ( v ac , ) of different 
transport means as v = v ac /(l + k), and expected value of the supply lead time 
as E(LT) = L/v; 

e. supply lead time ( LT) is a continuous random variable normally distributed. 

The logistic cost function considered in case of stochastic variability of supply 
lead time is: 



0 L = 0 () + 0 H + 0 T + 0 S + 0 EX [€/year] 



(3.40) 
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As already discussed in Sect. 2, in case of stochastic variability of the supply 
lead time ( LT ), three different situations may occur: 

(D) LT<E(LT) 

(E) E{LT) <LT<LT* 

(F) LT> LT* 

with LT* the maximum value of the supply lead time not causing a stock out 
event at a given service level ( SL ), and 



SL = prob(p TOT < LT* ■ d)= LT \ pdf(LT)dLT = F(LT*) (3.41) 

-00 



As a consequence, in order to evaluate the holding costs in Eq. (3.40) the 
expected inventory level ( E(1 )) in one ordering cycle has to be computed. With 
E(I) a , E(I) b , and E(I) C the expected inventory level in the three cases (Table 3.26), 
the expected inventory level during one ordering cycle can be evaluated as: 



E(LT ) 

E(I) = J E(I) a ■ pdf(LT)dLT 

LT* +=° 

+ J E(I) b ■ pdf(LT)dLT + \E{I) C ■ pdf(LT)dLT ■ 

E(LT) LT* 



(3.42) 



Table 3.26 Expected inventory level and ordering cycle length in the three cases considered 



Expected inventory level E(I) 



Case A 
LT < E(LT) 

Case B 

E(LT)<LT<LT* 



E(I) A =^Q + D- [E(LT) -ltj + ss 
E(I) b =—Q-D- [LT -E(LT)] + SS 



Case C 
LT> LT* 


E(I) c = l -Q + ^SS- 


1 [LT*-E(LT)J 
2 Q CT c 




Consumption time CT 


Case probability 


Case A 
LT < E{LT) 


Q 

D 


\pdj\z)dz 

—00 


Case B 

E(LT)<LT<LT* 


Q 

D 


\pdf{z)dz 

0 


Case C 
LT>Lf 


0- + ILT-LT* ] 
D 


+00 

J pdf ( z)dz = \ — SL 

z* 
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Under the assumption of E(I) C = H-Q (see Sect. 2), the expected inventory 
level E(I) in one ordering cycle can be computed as: 

E(I)= X -Q+Dg lt ■ [pdf (z*) + z* -F(z*)\ , (3.43) 

and the corresponding holding cost is: 

=c H -\^Q + D-(J LT -\pdf(z*)+z*-F{z*)^. (3.44) 

In Eqs. (3.43) and (3.44), z is the standardized variable of the stochastic 
variable LT , and hence: 

-* = LT * -E(LT') (3.45) 

a LT 

As far as concern shortage costs, they occur only in case C: 

= c s D ■ o it ' L(z*) , (3.46) 

with 

L(z*) = \z-z*)-pdf(z)dz= pdf(z*)-z*[l-F(z*)]. (3.47) 

z* 

In the model, quadratic functions are adopted in order to shape the dependency 
of both the transport cost and the external cost of transport on loss factor values: 

0 T (f,L) =Gm\af 2 +bf + c], ( 3 . 48 ) 

0 EX =GmL\Sf 2 +//•/]. (3.49) 

The logistics cost function can be finally reformulated as: 
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(3.50) 



and the corresponding logistic costs factor as 




\pdf{z*) + z*-F{z*)\^ + 



+ — ■ — + — • m ■ \ a ■ f 2 +b - f + cl 

c k Q c h L 7 7 J 




C H Q 



(3.51) 



For an assigned set of ( SL , cv LT , L) values, with a LT = cv LT ■ E(LT ), solving 
problem (3.29) the optimal values of loss factor (f OP i), order quantity (SOQ), 
reorder level ( r(f OPl )), and of the safety stock (SS(f OPT )) are obtained. 

The following procedure has been adopted here to solve problem (3.26): 

i. compute LT* from (3.41); 

ii. compute z* from (3.45); 

iii. compute f OP T and Popt by solving (3.29); 

iv. compute E(LT) = L/\>(f OPT y, 

v. compute r(fo P r) = D • E(LT); 

vi. compute SOQ = r(f OPT ) 'Popt, 

vii. compute a LT = cv LT ■ E(LT); 

viii. compute SStfopr) by means of: 



Step (iii) has to be carried out by means of a numerical method. 

Case Study 

The model has been applied to the case study discussed in Sect. 3.2.1 [9, 10]. 
Following the procedure above, the logistics problem (3.29) has been solved for 
PX product in site 1 (G = 9224 [unit/year], H = 3520 [h/year], D = 2.62 [unit/h]). 
Cost data in [€2013] have been considered ( c H = 18.98 [€ 2 oi 3 /unit-year], c s = 12.32 
[€ 2 oi 3 /unit]). A unit order cost ( c 0 ) of 100 [€/order] has been assumed. The 
corresponding EOQ value is 312 [unit/order]. 

Loss factor and average speed of transport values as in Sect. 3.2.1 and unit 
external cost data of water, rail, road, and air freight transport as in Sect. 3.1.2 
([€ 2013 / 0 km]) have been adopted. 

Parameter values adopted to evaluate unitary transport cost are in Table 3.27. 



SS = D-z* ■ cr LT 



(3.52) 
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Table 3.27 Parameters values for unitary transport cost evaluation per transport distance 



L 


a 


b 


c 


km 


re/kg] 


re/kg] 


[e/kg] 


200 


391.37 


-402.35 


108.69 


500 


329.84 


-259.17 


107.63 


1000 


227.28 


-20.54 


105.87 



EOQ and SOQ Comparison (In Case of Deterministic Supply Lead Time ) 

The hypothesis of deterministic supply lead time (cv = 0) allows to perform a 
preliminary comparison between results of the SOQ model proposed with those 
obtained by the traditional EOQ model [15]. 

The results obtained without considering cost of externalities (“ econ .”) are 
compared with results obtained considering external costs (“yz«7.”) (Table 3.28). 

In the sustainable solution, a SOQ value greater than the EOQ value is 
identified as optimal choice for each transport distance considered. When the 
transport distance considered increases, the optimal loss factor values decrease 
both in case of economic and in case of sustainable solutions. As a consequence, 
in case of long transport distance higher reorder levels are required. For the 
sustainable solution, when compared with the economic one, slower transport 
means (and higher reorder levels) are identified as optimal choice. 

Finally, as expected, sustainable solutions are characterized by higher logistics 
cost function (F L ) values than the corresponding ones in the economic case. 



Table 3.28 Sustainable and economic solutions comparison in case of cv = 0 



L 

km 


f sust. 


SOQ sust. 
unit 


Fl sust. 


r sust. 
unit 


200 


0.149 


318 


0.071 


12 


500 


0.037 


318 


0.073 


27 


1000 


0.037 


330 


0.074 


55 


L 

km 


f econ. 


SOQ econ. 
unit 


Fl econ. 


r econ. 
unit 


200 


0.335 


308 


0.069 


8 


500 


0.335 


308 


0.071 


20 


1000 


0.037 


318 


0.073 


55 



Evaluating the Effects of the Supply Lead Time Variability 
In case of a variable supply lead time, shortage costs affect solution of (3.29). 
Shortage costs depend on safety stock level ( SS ), which is a function of the lead 
time variability and of the SL (see Eq. 3.52). In order to investigate the magnitude 
of the effects of the supply lead time variability on the solutions of the logistics 
problem (3.29), it has been solved for different values of the lead time variability 
(cv) and of the service level (SL). A negligible influence of the service level values 
on the solution of the logistics problem (Jopt, SOQ, r(f 0 pT)) has been observed. 
This is due to the small value of the unitary shortage cost ( c s ) considered. 
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As an example, results obtained in case of SL = 0.95 are in Table 3.29 and 
Table 3.30. Road transport is identified as optimal choice for each level of supply 
lead time variability considered in case of short transport distance (L = 200 [km]). 
Fast road transport means are solution of problem (3.29) when a high variability of 
the lead time is considered. They are solution of (3.29) for higher supply lead time 
variability value in case of long transport distance (L = 500-1000 [km]). 
Transport means characterized by better environmental performances (rail, ship) 
are the optimal choice in the other cases considered (see Table 3.29). Reorder 
level obtained ( r(f OPT ), see Table 3.30) are consistent with the loss factor values 
identified as optimal by means of (3.29): smaller reorder level can be adopted in 
case of faster transport means. In case of short transport distance and low 
variability of supply lead time, SOQ values higher than EOQ are observed. 
Furthermore, supply lead time variability weakly affects SOQ values observed. 
On the contrary, an increase in the supply lead time variability (and in the 
transport distance) leads to higher SS values (see Table 3.30). The effect of the 
supply lead time variability on total logistic costs ( F L (f 0 pT :), see Table 3.30) 
proved to be dependent on the transport distance considered. In the two limit cases 
considered (cv = 0.00 and cv = 2.00, SL = 0.95) an increase in total annual costs 
of 5.5 %, 18.7 %, and 45.7 % is observed in case of a transport distance of 200 
[km], 500 [km], and 1000 [km], respectively (see Table 3.30). 



Table 3.29 Optimal means of transport (/opt) and SOQ values for different L and cv values in 

case of SL = 0.95 







foPT 








SOQ 




CV 


200 


L [km] 
500 


1000 


CV 


200 


L [km] 
500 


1000 


0.00 


Truck 7.5-12 [t] 


Rail — electric 


Rail — electric 


0.00 


318 


318 


330 


0.10 


Truck 7.5-12 [t] 


Rail — electric 


Ship 


0.10 


318 


318 


318 


0.25 


Vans 


Rail — electric 


Ship 


0.25 


319 


318 


318 


0.50 


Vans 


Rail — electric 


Ship 


0.50 


319 


318 


318 


0.75 


Vans 


Vans 


Ship 


0.75 


319 


319 


318 


1.00 


Vans 


Vans 


Ship 


1.00 


319 


319 


318 


2.00 


Vans 


Vans 


Vans 


2.00 


319 


319 


318 
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Table 3.30 Reorder level (r(fopr)), SS, and Fp values for different L and cv values in case of 

SL = 0.95 



r (fopr) SSffopr) 



CV 


200 


L [km] 

500 1000 


CV 


200 


L [km] 

500 1000 


0.00 


25 


55 


109 


0.00 


0 


0 


0 


0.10 


25 


55 


102 


0.10 


2 


4 


8 


0.25 


7 


55 


102 


0.25 


2 


11 


21 


0.50 


7 


55 


102 


0.50 


3 


22 


42 


0.75 


7 


19 


102 


0.75 


5 


12 


63 


1.00 


7 


19 


102 


1.00 


6 


15 


84 


2.00 


7 


19 


37 


2.00 


12 


31 


61 



F L (foPT Popt) 



CV 


200 


L [km] 
500 


1000 


0.00 


0.070 


0.073 


0.074 


0.10 


0.071 


0.074 


0.076 


0.25 


0.071 


0.076 


0.080 


0.50 


0.072 


0.079 


0.085 


0.75 


0.072 


0.081 


0.091 


1.00 


0.073 


0.082 


0.097 


2.00 


0.074 


0.087 


0.107 



In order to investigate the effect of the external costs on the solution of 
problem (3.29), it has been solved in case no external costs are computed in the 
logistics cost function (@ex = 0). Results obtained in case of a service level (SL) 
equal to 0.95 are in Table 3.31. 

Table 3.31 Optimal loss factor (/opt ) values of the sustainable and of the economic 

( ( Lex = 0) solution 



fppT SUSt 



CV 


200 


L [km] 
500 


1000 


0.00 


0.149 


0.037 


0.037 


0.10 


0.149 


0.037 


0.010 


0.25 


0.554 


0.037 


0.010 


0.50 


0.554 


0.037 


0.010 


0.75 


0.554 


0.554 


0.010 


1.00 


0.554 


0.554 


0.010 


2.00 


0.554 


0.554 


0.554 



foPT econ. 



CV 


200 


L [km] 
500 


1000 


0.00 


0.335 


0.335 


0.037 


0.10 


0.335 


0.335 


0.010 


0.25 


0.554 


0.554 


0.554 


0.50 


0.554 


0.554 


0.554 


0.75 


0.554 


0.554 


0.554 


1.00 


0.554 


0.554 


0.554 


2.00 


0.554 


0.554 


0.984 
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As showed in Table 3.31, the internalization of the external costs leads to 
identify as optimal choice slower transport means (3.29). As far as concern SOQ 
values, negligible differences have been observed in the solutions obtained in the 
two cases. 

A sensitive analysis has been carried in order to evaluate the effects of a 
change in the ordering cost on the solutions of (3.29). For this purpose, problem 
(3.29) has been solved considering other two unitary order cost values (10 
[€/order] and 200 [€/order]). Results obtained are in Table 3.32. Solutions 
obtained proved to be slightly affected by order cost values, except for SOQ 
values: with the increase of fixed order cost, higher values of SOQ minimize total 
logistic costs (see Table 3.32 and Table 3.29). 



Table 3.32 SOQ values for different L, cv, and Co values in case of SL = 0.95 



SOQ 


Co = 


= 10 [e/order] 


Co = 


200 [€/order] 


CV 


L [km] 
200 


500 


1000 


L [km] 
200 


500 


1000 


0.00 


93 


98 


100 


403 


420 


420 


0.10 


93 


99 


100 


403 


420 


420 


0.25 


93 


99 


98 


403 


420 


464 


0.50 


93 


98 


109 


374 


420 


464 


0.75 


93 


98 


109 


374 


467 


464 


1.00 


93 


104 


104 


374 


467 


464 


2.00 


93 


104 


98 


374 


467 


467 



